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Spacetime Algebra (STA) provides unified, matrix-free spinor methods for rotational dynamics 
in classical theory as well as quantum mechanics. That makes it an ideal tool for studying particle 
models of zitterbewegung and using them to study zitterbewegung in the Dirac theory. This paper 
develops a self-contained dynamical model of the electron as a lightlike particle with helical zitter- 
bewegung and electromagnetic interactions. It attributes to the electron an electric dipole moment 
oscillating with ultrahigh frequency, and the possibility of observing this directly as a resonance in 
electron channeling is analyzed in detail. A modification of the Dirac equation is suggested to incor- 
porate the oscillating dipole moment. That enables extension of the Dirac equation to incorporate 
electroweak interactions in a new way. 

Keywords: zitterbewegung, geometric algebra, electron channeling, de Broglie frequency 



I. INTRODUCTION 

This paper continues a research program investigating 
implications of the Real Dirac Equation for the interpre- 
tation and extension of quantum mechanics. Details of 
the program have been reviewed elsewhere [1-3], so it 
suffices here to state the main ideas and conclusions to 
set the stage for the present study. 

The program began with a reformulation of the Dirac 
equation in terms of Spacetime Algebra (Section II), 
which revealed geometric structure that is suppressed in 
the standard matrix version. In particular, it revealed 
that the generator of phase and electromagnetic gauge 
transformations is a spacelike bivector specified by elec- 
tron spin. In other words, spin and phase are inseparably 
related — spin is not simply an add-on, but an essential 
feature of quantum mechanics. However, physical impli- 
cations of this fact depend critically on relations of the 
Dirac wave function to physical observables, which are 
not specified by the Dirac equation itself. That started 
the present research program to investigate various pos- 
sibilities. 

A standard observable in Dirac theory is the Dirac cur- 
rent, which doubles as a probability current and a charge 
current. However, this does not account for the mag- 
netic moment of the electron, which many investigators 
conjecture is due to a circulation of charge. But what is 
the nature of this circulation? After a lengthy analysis 
of the Dirac equation Bohm and Hiley conclude [4] : "the 
electron must still be regarded as a simple point particle 
whose only intrinsic property is its position." Under this 
assumption, spin and phase must be kinematical features 
of electron motion. The charge circulation that generates 
the magnetic moment can then be identified with the zit- 
terbewegung of Schroedinger [5]. 

This raises the central question of the present research: 
Is the zitterbewegung, so construed, a real physical phc- 
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nomenon, or is it merely a colorful metaphor? Although 
this question was motivated by structural features of the 
Dirac equation, it cannot be answered without attribut- 
ing substructure to electron motion that is not specified 
by standard Dirac theory. 

The main purpose of this paper is to formulate and 
study a well-defined particle model of the electron with 
spin and zitterbewegung dynamics motivated by the Dirac 
equation. Since the term zitterbewegung is quite a mouth- 
ful, I often abbreviate it to zitter, especially when it is 
used as an adjective. 

We study the structure of the zitter model in consider- 
able detail with the aim of identifying new experimental 
implications. The main conclusion is that the electron 
is the seat of a rapidly rotating electric dipole moment 
fluctuating with the zitter frequency of Schroedinger. As 
this frequency is so rapid, it is observable only under res- 
onance conditions. It is argued that many familiar quan- 
tum mechanical effects may be attributable to zitter res- 
onance. Moreover, the new possibility of observing zitter 
directly as a resonance in electron channeling is analyzed 
in detail, because the prospects of crucial experimental 
tests are very promising. 

The relation of the zitter particle model to the Dirac 
equation is also studied. The main conclusion is that, 
though zitter oscillations are inherent in the Dirac equa- 
tion, they will not be manifested as an oscillating electric 
dipole without altering the definition of charge current. 
A simple modification of the Dirac equation to incorpo- 
rate the altered definition is proposed. Remarkably, that 
opens the door for incorporating electroweak interactions 
into the Dirac equation in a novel way. 

In conclusion, the relation of the zitter particle model 
to the Dirac equation can be considered from two differ- 
ent perspectives. On the one hand, it can be regarded 
as a "quasiclassical" approximation that embodies struc- 
tural features of the Dirac equation in a convenient form 
for analysis. On the other hand, it can be regarded as for- 
mulating fundamental properties of the electron that are 
manifested in the Dirac equation in some kind of average 
form. The choice of perspective is left to the reader. 
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II. SPACETIME ALGEBRA 

Spacetime algebra is thoroughly expounded elsewhere 
[1], so a brief description is sufficient here, mainly to de- 
fine terms. 

We represent Minkowski spacetime as a real 4- 
dimensional vector space M. 4 . The two properties of 
scalar multiplication and vector addition in A4 provide 
only a partial specification of spacetime geometry. To 
complete the specification we introduce an associative ge- 
ometric product among vectors a, b, c, . . . with the prop- 
erty that the square of any vector is a (real) scalar. Thus 
for any vector a we can write 



aa 



(1) 



where e is the signature of a and \a\ is a (real) positive 
scalar. As usual, we say that a is timelike, lightlike or 
spacelike if its signature is positive (e = 1), null (e = 0), 
or negative (e = — 1). We can specify the signature of A4 4 
as a whole, by adopting the axioms: (a) M 4 contains at 
least one timelike vector; and (b) every 2-plane in M 4 
contains at least one spacelike vector. 

To facilitate applications of STA to physics a few def- 
initions and theorems arc needed. From the geometric 
product uv of two vectors it is convenient to define two 
other products. The inner product u • v is defined by 

u • v = | (uv + vu) = v • u , (2) 
while the outer product u A v is defined by 

u Av — \(uv — vu) = — v A u . (3) 
The three products are therefore related by 



uv = u • v + u A v . 



(4) 



This can be regarded as a decomposition of the product 
uv into symmetric and skewsymmetric parts, or alterna- 
tively, into scalar and bivector parts. 

The inner and outer products can be generalized. One 
way is to define the outer product along with the notion 
of k-vector iteratively as follows: scalars are defined to 
be 0-vectors, vectors are 1-vectors, and bivectors, such as 
u A v, are 2-vectors. For a given fc-vector K, the integer 
k is called the grade of K. For k > 1, the outer product 
of a vector v with a k- vector A" is a (k + 1 )-vector defined 
in terms of the geometric product by 

v A K = \{vK + (—l) k Kv) = (—l) k K A v . (5) 

The corresponding inner product is defined by 

v K = \(vK + (-l) k+1 Kv) = (-l) k+1 K -v, (6) 

and it can be proved that the result is a (k — f )-vector. 
Adding (5) and (6) we obtain 



vK = v • K + v A K , 



(7) 



which obviously generalizes (4). The important thing 
about (7), is that it decomposes vK into (k — l)-vector 
and (k + 1 )-vector parts. 



By continuing as above, STA as been developed into 
a complete coordinate- free calculus for spacetime physics 
[1]. However, to hasten comparison with standard Dirac 
algebra, we interrupt that approach to introduce coor- 
dinates and a basis for the algebra. Let {7^; 0,1,2,3} 
be a right-handed orthonormal frame of vectors with 70 
in the forward light cone. In accordance with (2), the 
components g of the metric tensor for this frame are 
given by 



S/^!/ In ' 7^ — 2 (t^TV + 7^7^) • 



(8) 



This will be recognized as isomorphic to a famous formula 
of Dirac's. The difference here is that the 7^ are vectors 
rather than matrices. A coframe {7^} is defined by the 
usual convention for raising and lowering indices: 7^ = 

The unit pseudoscalar i for spacetime is related to the 
frame {7^} by the equation 



i = 7o7i7273 = 7o A 71 A 72 A 73 . 



(9) 



It is readily verified from (9) that i 2 = —1, and the geo- 
metric product of i with any vector is anticommutative. 

By multiplication the 7^ generate a complete basis of 
/c-vectors for STA, consisting of the 2 4 = 16 linearly in- 
dependent elements 



1, 7m> 7m a lv, 7nh i- 



(10) 



Multivectors with even grade compose a subalgebra of 
the STA generated by the bivectors {cr fc = 7^70; & = 
1,2,3}, so that 



ci£r2<x 3 = 70717273 = 



(11) 



Any multivector can be expressed as a linear combination 
of these elements. 

For example, a bivector F has the expansion 



F=\F^ ltl A lvi 



(12) 



with its "scalar components" P" given by 

F» v =Y-F-Y =Y-{^-F) = (7^7^-^.(13) 

Note that the two inner products in the middle term 
can be performed in either order, so a parenthesis is not 
needed; also, we use the usual convention for raising and 
lowering indices. 

Alternatively, one can decompose the bivector into 



F = E + zB = F ko cr k + \F k ^CTj(T k , 



(14) 



corresponding to the split of an electromagnetic field 
into "electric and magnetic parts." 

The entire spacetime algebra is obtained from linear 
combinations of basis k- vectors in (10). A generic ele- 
ment M of the STA, called a multivector, can thus be 
written in the expanded form 



M = a + a + F + bi + (3i = ^(M) k , 



(15) 



k=0 
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where a and (3 are scalars, a and b are vectors, and F 
is a bivector. This is a decomposition of M into its li- 
vedo? parts, with k = 0, 1, 2, 3, 4, where (. . means "fc- 
vector part." Of course, (M) = a, (M)i=a, (M) 2 =F, 
(M) 3 = W, (M} 4 = 

Computations are facilitated by the operation of rever- 
sion. For M in the expanded form (15) the reverse M 
can be defined by 

M = a + a-F-bi + (3i. (16) 

Note, in particular, the effect of reversion on the various 
k- vector parts: 

a = a, a = a, F = —F, i = i . It is not difficult to 
prove that 

(MN)~=NM, (17) 

for arbitrary multivectors M and N. For scalar parts 
(M) = (M) , we have 

(M) = (M) hence (MN) = (NM) . (18) 

The scalar part in STA corresponds to the trace in Dirac's 
matrix representation of the algebra. 

Besides the inner and outer products defined above, 
many other products can be defined in terms of the ge- 
ometric product. We will need the commutator product 
M x N, defined for any M and N by 

M x N = \{MN - NM) — —N x M . (19) 

We will make frequent use of the "Jacobi identity" in the 
form of a "derivation" : 

M x{N x P) = (M x N)x P + N x{M x P) .(20) 

For bivectors S and F we can now expand the geometric 
product as follows: 

SF = S-F + SxF + SAF (21) 

where S x F is a bivector and S A F is a pseudoscalar. 
Note how this differs from the expansion (7). 

One great advantage of STA is that it provides 
coordinate-free representations of both tensors and 
spinors in the same system. Spinors can be represented as 
even multivectors serving as algebraic operators. Using 
reversion, it is easy to prove that every even multivector 
tp satisfies 

W = pe if) . (22) 
Hence, for p ^ 0, ip can be written in the canonical form 

V> = (f>e if> )*R, (23) 

where 

RR = RR = 1 . (24) 

Every such R determines a Lorentz rotation of a given 
multivector M: 

R : M -> M' = RMR , (25) 



and every Lorentz rotation can be expressed in this 
coordinate-free form. Construed as an operator in this 
sense, the quantity R is called a rotor while ip is called 
a spinor. We shall describe particle kinematics and the 
Dirac wave function by spinors in this sense. 

The set of all rotors form a multiplicative group called 
the rotor group or the spin group of spacetime. When 
R = R(t) represents a one parameter family of rotors 
(hence of Lorentz rotations), "angular velocities" O and 
fy are defined by 

R = ±SIR = \RQ! , (26) 

or 

fl = 2RR = R fl'R and fl' = 2RR , (27) 

where the overdot indicates derivative. It follows from 
(24) and (15) that and 17' are bivectors. 

We represent each spacetime point as a vector x — 
x^ln with rectangular coordinates = 7^ • x. The 
vector derivative with respect to the point x is defined 
by V = d x = -f^dxn. As V is a vectorial operator, we 
can use (7) to decompose the derivative of a fc-vector field 
K = K(x) into 

VK = V ■ K + V AK , (28) 

where the terms on the right can be identified, respec- 
tively, with the usual divergence and curl in tensor cal- 
culus. 

Besides the STA definitions and relations given above, 
many others can be found in the references. We shall 
introduce some of them as needed. 



III. CLASSICAL PARTICLES WITH SPIN 

Classical models of the electron as a point particle with 
spin were first formulated by Frenkel [6] and Thomas [7] , 
improved by Mathisson [8] , and given a definitive form by 
Weyssenhoff [9]. They have been revisited from time to 
time by many investigators [10-12], including others to 
be cited below. They are of interest mainly for the insight 
they bring to the interpretation of quantum mechanics. 
But they are also of practical value, for example, in the 
study of spin precession [1, 13] and tunneling [14-16]. 

In Weyssenhoff 's analysis the models fall into distinct 
classes, differentiated by the assumption that the elec- 
tron's spacetime history is timelike in one and lightlikc 
in the other. The timelike case has been studied by many 
investigators, such as those cited above. Ironically, the 
lightlike case, which Weyssenhoff regarded as far more 
interesting, has been generally ignored. Without being 
aware of his analysis, I arrived at similar conclusions 
from my study of real Dirac theory. This paper revis- 
its Weyssenhoff theory with new mathematical tools to 
simplify and extend it. 

Our first task is to reformulate Weyssenhoff theory in 
the language of spacetime algebra (STA) and show how it 
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is simplified, clarified and extended. A decided advantage 
of STA is its uniform use of spinors in both classical and 
quantum theory. For one thing, this makes it easier to 
relate classical models to solutions of the Dirac equation 
[1]. Although we are most interested in the lightlike case, 
we treat the timelikc case in enough detail to compare 
the two cases. Besides, the timelike models can be used 
for other particles with spin besides electrons, including 
atoms. 

We consider a particle with spacetime history z = z(t) 
and proper velocity [1] 



dz 
J?' 



(29) 



For the time being, we allow u to be cither timelike, with 



1, or lightlike, with u — 0, so in either case 



u • u = . 



(30) 



In the timelike case, of course, the parameter r is taken 
to be proper time, but in the lightlike case it must be 
regarded as arbitrary for the time being. 

Suppose the particle has proper momentum p = p(r) 
and intrinsic angular momentum (spin) represented by 
a bivector S = S(t). Equations of motion are obtained 
from general conservation laws for momentum and an- 
gular momentum. However, to formulate these laws cor- 
rectly it is essential to note that just as p contributes 
an orbital momentum, S contributes an intrinsic part to 
the momentum, so we cannot make the usual assumption 
p = mu. Instead, the relation between p and u depends 
on the dynamics of motion. The same can be said about 
the relation between S and u, although we shall assume 
that it is restricted by the kinematical constraint 



S-u = 0. 



(31) 



This reduces the degrees of freedom in S to three at most. 

The noncollinearity of p and u raises a question about 
how mass should be defined. Without prejudicing the 
issue, it is convenient to introduce a dynamical mass de- 
fined by 



m=p-u. 



(32) 



This quantity is well defined for both timelike and light- 
like particles, so lightlike particles can have mass in this 
sense. However, the value of this mass obviously depends 
on the choice of time parameter r, which, so far, is not 
well defined in the lightlike case. This is our first hint 
that mass is intimately related to intrinsic time scaling 
of electron histories. 

Without presupposing a relation between momentum 
and velocity, momentum conservation can be given the 
general form 



P = f, 



(33) 



where the proper force / describes transfer of momentum 
through external interactions. 



Angular momentum conservation is governed by 

S = uAp + T, (34) 

where the proper torque T describes angular momentum 
transfer through external interactions. To see where (34) 
comes from, we introduce the orbital angular momentum 
pAzso the total angular momentum is 



J = pAz + S . 



(35) 



The angular momentum conservation then has the more 
familiar form 



j = f a z + r . 



(36) 



where f Az is the orbital torque. The equivalence of (36) 
with (34) follows immediately by differentiating (35). 

Equation (36) tells us that J is a constant of motion for 
a free particle with spin. Spin is not separately conserved. 
According to (35), angular momentum can be exchanged 
back and forth between orbital and spin parts. 

To get well defined equations of motion from the con- 
servation laws we need to specify the interactions. For 
a particle with charge q and dipolc moment M = M(t), 
we use 



/ = qF -u + VF • M 



(37) 



where F = F(z) is an applied electromagnetic (bivector) 
field [aptly called the Faraday] and the vector derivative 
V = 7^9^ operates on F. The first term on the right of 
(37) is the standard "Lorentz force," while the second is a 
force of "Stern-Gerlach" type. The latter is accompanied 
by the torque 



F x M . 



(38) 



The "interaction laws" (37) and (38) must be supple- 
mented by a "constitutive equation" expressing M as a 
function of S and u. In general M can include both 
electric and magnetic dipole moments. The general con- 
stitutive constraint 



M x S = 



(39) 



is satisfied by the standard constitutive equation M = kS 
for a magnetic moment, but also by the more general 
relation 



M = nSe i& . 



(40) 



The factor e l/3 is a "duality rotation" relating a mag- 
netic moment iSn sin ip to an electric moment Sn cos ip. 
As noted later, the magnetic moment (density) in Dirac 
theory has this form. In the lightlike case we shall see 
that there is an alternative way to satisfy (39). 

From the specific form of the interaction laws we can 
prove that the magnitude of the spin is a conserved quan- 
tity. First note that (31) implies 



SAS = 



so 



S-S. 



(41) 
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This follows from the identity 

(S A S)u = (S A S) • u = 2S A (S • u) . 
From (34) with (38) we get 

i A-S 2 = S ■ S = S • (u Ap) + S ■ (M x F) . 

The first term on the right vanishes because 

S ■ (u A p) = (5 • u) ■ p = , 
and the second term vanishes by (39) and the identity 

S ■ (F x M) = (SFM) =F-(M x S) . (42) 
Therefore 



S 2 = - \S\ 2 = constant. 



(43) 



The negative sign appears because (31) implies that S 
cannot be a timelike bivector. 

Up to this point, our equations apply equally to time- 
like and lightlike particles. To go further we must con- 
sider each of these cases separately. There are three dis- 
tinct possibilities which have been studied in the past. 
Most investigators have been attracted by the possibility 
of choosing \S\ = Ti/2 to make contact with Dirac theory. 
The possibility of coupling this with a timelike velocity 
has been most thoroughly studied by Corben [10]. How- 
ever, his model does not recognize the intimate connec- 
tion between electron spin and phase. The alternative 
possibility of coupling with a lightlike velocity overcomes 
that deficiency. I confess to struggling with that model 
for many years before I realized that for a lightlike par- 
ticle the spin must be a lightlike bivector. Shortly there- 
after I learned that Weyssenhoff had figured that out long 
before [9]. 

There are many more things to say about the timelikc 
case, but we shall concentrate on the more interesting 
lightlike case of Weyssenhoff. The tools of STA enable us 
to develop this case much farther than he could, indeed, 
to a fairly definitive and satisfactory conclusion. 



IV. FREE PARTICLE MOTION 

We gain crucial insight into the "kinematics of spin" 
from the free particle solution. The solution has the same 
form for both timelike and lightlike cases, though we shall 
see that there are significant differences in physical inter- 
pretation. 

For a free particle it follows from (33), (36) and (35) 
that p and J = pAz + S are constants of motion. We seek 
a general solution with constant mass m e = p • u = p • z. 
Integrating this relation, we obtain 



p • (z - zq) = m e T , 



(44) 



where zo is the particle position at r = 0. Angular mo- 
mentum conservation gives us 



where So is the initial value of the spin. Adding these 
equations and multiplying by p^ 1 = p/p 2 we get an equa- 
tion for the particle history: 

z = [S(T)-S ]p- 1 +m e p- 1 T + z . (46) 
We can write this in the form 

z(t) = r(r) + x{t) , (47) 
where a center of motion is defined by 

x = Utl T + z - So • p- 1 , (48) 
P 

and a radius vector is defined by 

r = S(T)-p- 1 . (49) 
Differentiating, we obtain 

r = S • p^ 1 = (u A p) • p^ 1 = u — niep^ 1 . (50) 

It follows that r*r = 0, so r 2 is a constant of motion, and 
r = r(r) must be a rotating vector. Hence (47) describes 
a helix in spacetime centered on the timelikc straight line 
(48). 

We get an explicit form for the temporal behavior of 
S and r by introducing a rotor R = R(t) so that 



S — RSqR . 
Then (49) is consistent with 

r = Rr R 
provided 

p- Q = 
for f2 = 2RR. Furthermore 

z = u = RuqR 
if SI is a constant, in which case R' integrates to 



R 



e 2 



Sir 



(51) 
(52) 
(53) 
(54) 
(55) 



From this we conclude that all aspects of the motion are 
determined by the scalar, vector and bivector constants 
of motion, m e , p and 0. The values of these constants 
are left to be determined by other considerations. 

These results apply to both timelike and lightlike cases. 
In the timelike case, if one assumes |5| = U/2, it is easy 
to prove that 



mi 



2m e c 2 
h 



1.5527 x 10 21 rad/sec, (56) 



which is the zitter frequency of Schrocdinger where m e is 
the electron rest mass. 

Moreover, writing p = m e v, we can express the proper 
time in the form r = v • z = v • x, so (55) can be expressed 
in the form 



(z- z ) Ap = S(t) - S , 



(45) 



R = eM^ T = e ip - x/h . 



(57) 
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which, with i as a unit spacelike bivector, has been shown 
elsewhere [1] to correspond precisely to plane wave solu- 
tions of the Dirac equation. 

The same solution applies to the lightlike case, but 
with the additional constraint uj e X e — c, so the electron 
zitter radius has the fixed value 



A, 



n 



2m e c 



1.9308 x 1(T 3 A. 



(58) 



Regarded as a solution of the Dirac equation, the spinor 
(57) describes a congruence of either timelike or lightlike 
histories, depending on the choice of velocity observable. 
Our next task is to generalize these insights to include 
interactions. 



V. LAGRANGIAN FOR A LIGHTLIKE 
PARTICLE WITH SPIN 

In Section III we postulated equations of motion and 
a number of relations among velocity, spin, momentum 
and mass variables. To ensure that all this constitutes a 
complete, coherent and self-consistent dynamical system 
we show that it can be derived from a single Lagrangian, 
but with introduction of further constraints. 

Conventional Lagrangians cannot be used for lightlikc 
particles, owing to the absence of proper time as a nat- 
ural parameter. Weysscnhoff [17] addressed this prob- 
lem by deriving parameter invariant Lagrangians for both 
timelikc and lightlikc particles. However, he encountered 
problems that limited his treatment to the free particle 
case. Kriiger [18] resolved the problem for a lightlikc 
particle by introducing a Lagrangian defined in terms of 
particle path curvature (i.e. acceleration) instead of ve- 
locity. However, his model does not have all the features 
that we are looking for, so we employ a different approach 
here. 

We take full advantage of STA by using it to construct 
a spinor-based Lagrangian. Proca [19] was the first to use 
Dirac spinors for describing classical particles. Barut and 
Zanghi [20] used the same approach to model electron 
zitter. Gull [21] translated their model into STA and 
noted that it lacked a magnetic moment. 

Here we see that spinors are especially helpful for build- 
ing lightlike constraints into the Lagrangian. The result- 
ing spinor equations provide a superior dynamical model 
for a lightlike particle with spin, mass and zitter, which, 
as we see later, is intimately related to the Dirac equa- 
tion. Thus, we are able to complete Weyssenhoff 's pro- 
gram to construct dynamical equations for a lightlike par- 
ticle with zitter and electromagnetic interactions. 

Continuous "motion" of a particle in spacetime is rep- 
resented by a curve z = z(ip) and its derivative z' — 
dz/dip, where (p is an affine parameter for which the phys- 
ical interpretation is initially unspecified. The kinematic 
structure of this curve is described by a spinor ip — ip(ip) 
and its derivative ip'. 



The dynamics of motion in an external electromagnetic 
field F(z) — V A A with vector potential A = A(z) is 
determined by a Lagrangian L = L(z, ip, P) of the form 



l = ( -^'7+71^ + p( z ' - V^y+VO 

+qAip^ + ip + gA e F-07+7iV' ) > (59) 



with 7± = 70 ± 72 and units: c = 1,U = m c = fixed 
reference mass. The coupling constants are charge q and 
a length A e , which we anticipate identifying as the free 
electron zitter radius (58), so it amounts to introducing 
the electron rest mass m e as a coupling constant. A 
vectorial Lagrange multiplier P relates the two kinds of 
kinematic variable. 

The method of multivector differentiation is the sim- 
plest and most elegant way to derive equations of motion 
[22]. For a Lagrangian that is homogeneous of degree 
one in derivatives, variation with respect to a multivec- 
tor variable X = X((p) yields the multivector Lagrange 
equation 



5xL = d x L-d ip {dx>L) = 0, 



(60) 



where X' = d v X. 



Variation of the Lagrangian with respect to P obvi- 
ously gives 



z 1 = tp-y+ip = w , 



(61) 



which defines a particle velocity w in terms of the spinor 
The result of varying position vector z is the force law 



p' = qF -w + VF ■ M u 



(62) 



with momentum vector p defined by the canonical ex- 
pression 

p = P - qA, (63) 
and electromagnetic moment bivector defined by 

M w = qAe-07+71^. (64) 

Of course, the vector derivative V in (62) operates only 
on F and not on M. 

The result of varying spinor ip in the Lagrangian is the 
dynamical spinor equation 

?i-V/7+7i = -W7+ + gA e FV>7+7i . (65) 
From this one can easily derive the equation of motion 

S' w = w Ap + F x M w (66) 
for a spin bivector defined by 



Sw = 1^7+71^- 
One can also derive the equation 



-w 



h 



v Sw + qKF ■ W, 



(67) 



(68) 
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which relates particle acceleration to spin and momen- 
tum. 

The Lagrangian has delivered a system of coupled 
equations for particle velocity, spin and momentum. 
Some simplifications are still needed to facilitate phys- 
ical interpretation and analysis. 



VI. OBSERVABLES AND DYNAMICAL 
STRUCTURE 

Consider implications of the canonical decomposition 

ip = (pe l P)2R. The rotor R = R((p) determines a one 
parameter family of Lorentz rotations 



(69) 



that transforms a fixed orthonormal frame of vectors into 
an intrinsic comoving frame following the particle. The 
comoving frame is coupled to the particle velocity by 

w = z' = tp-f+ip = pu , (70) 
where a new rescaled velocity has been defined by 

u = Rj + R = e + e 2 . (71) 

Thus, we identify p as a time scale parameter! 

The duality factor e l @ rotates multivectors into their 
duals, as shown by the equation 



e J/3 7M = 7m cos/3 + 27^ sin/3 = 7^" 



-iff 



(72) 



Note that anticommutivity of vectors with the pscudo- 
calar has eliminated the duality factor from (70). How- 
ever, the duality factor commutes with bivectors. So for 
the spin bivector we find 



S w = ^ 1+ll ^^^pR 1+ RR 1 iRe^ 
= \pue^). 



(73) 



This makes the relation of spin to velocity u = e + e 2 
explicit. Further simplification is possible by eliminating 
the duality factor. 

Note that the unit pseudoscalar can be written i — 
7o7i7273 = e eie 2 e 3 = e e 2 e 3 ei, so ue e 2 = u implies 
iu = —ui = ueie^. Consequently, the effect of the duality 
factor multiplying a null vector is equivalent to a rotation, 
and (73) is proportional to the form 



e %t3 ue\ 



ue~ i0 e x = ue ieoe2/3 ei = ue e3£l/3 ei 



= ?i(ei cos/3 + e 3 sin/3). 



(74) 



In other words, the net effect of e l/3 in (73) is to rotate 
ei into e 3 . Since the duality factor appears in every term 
of the spinor equation (65), it can be eliminated by ab- 
sorbing it into the rotor, as shown by 



■07-t 



) 2i?e^ /2 7+ = p2(i?e 7l73/3/2 ) 7+ , 



so that (73) reduces to 

S w = |"07+7i^ = = pS, 



(75) 



where we have introduced a rescaled spin bivector S. 
Note that the spin has a dual form 

-lb = -e 3 eiuei = -^e 3 u, 

Thus, we can characterize the spin bivector by the com- 
pact equations 

S = -^ ue i = tsu , where s = — e 3 (76) 
defines a spin vector. 



A. Momentum, mass and velocity 

Particle mass relates momentum to velocity. When 
momentum is not collinear with velocity there are two 
distinct concepts of mass. First, assuming that the ki- 
netic momentum p points to the future, we can define a 
positive mass m p in the usual way: p 2 = m 2 > 0. Sec- 
ond, projection of momentum p onto the particle velocity 
u defines a dynamical mass m = p • u. This turns out to 
be the most significant mass variable. 

To specify the relation of momentum to velocity more 
completely, we note that the comoving frame {e M } is not 
uniquely determined by its relation to the particle veloc- 
ity in (71) and spin in (76). Consequently, we are free to 
constrain e to lie in the p Au plane, as specified by the 
equation 



p A u A e = p A e 2 A e = 0. 



(77) 



Hence, we can introduce mass parameters mi and m 2 so 
that 



p = mie - m 2 e 2 , 

with m = p • u = mi + m 2 . 

The various masses are thus related by 

^2 ra 2 ™2 



mj — m 2 mi — m 2 



Hz - r 

m 2 (p • u) 2 (mi + m 2 ) 2 mi + m 2 



(78) 



(79) 



Functional values for the masses will be determined by 
the equations of motion. 

Now the relation of momentum to particle velocity can 
be expressed in the simple form 

up = m(l + e 2 e ) with uAp = me 2 eQ. (80) 

Likewise, it follows from (78) thatp-ei andp-e 3 = p-s = 
0. Hence, the relation of momentum to spin is reduced to 
the simple equations 



p * o — —p • (uei) = -met, 



(81) 



p A S = p A (isu) = ip ■ (us) = ims 
or ms = —i(p A S) = p ■ (iS). 



(82) 



These relations considerably simplify the velocity equa- 
tion of motion (68) and, as we see later, the spinor equa- 
tion (65) by reducing direct coupling with momentum to 
the scalar variable m = p • u. One more simplification is 
needed. 



B. Time scaling and gauge transformations 

For purposes of measurement, we must relate the affinc 
parameter ip to some observable time scale. This can 
be done by projecting the electron's null velocity onto a 
timelike reference history x — x(t) with proper time r 
and velocity 



v = x = dx/dr, 



(83) 



thus associating a proper time with the electron motion. 
The electron velocity (6f ) is thus rescaled to proper time 
by 

dr dr dz 

p = v-w = —, so w = - — — = pu. (84) 
dip dip dr 

Accordingly, the rescaled particle velocity is given by 



dz 
dr~ 



(85) 



This leaves us with the problem of specifying the ref- 
erence history and determining the time scale factor p. 
We have considerable leeway in choosing the reference 
history, the only requirement being that its relation to 
the electron history z(r) be well defined. Two obvious 
choices for reference velocity are v = p/m p or u = eo, 
but we shall see that they have drawbacks, so it is best 
to leave the option open until a decision is needed. 

Recall that ip was originally introduced as an arbitrary 
affine parameter in the Lagrangian (59). In fact, the 
Lagrangian time scale depends on the spinor gauge, as 
we now prove. Dropping the prime notation for d/dip 
from here on, and using the prime to indicate a change of 
variable, we note that a time scale gauge transformation 
of the spinor can be defined as a boost 



which induces a scale change 



(86) 



-07 + V ~~ > ip'l+ip' = e Q, 07+'0- (87) 
Combining this with the change of variable 

d d dip d , . 



we have a time scale gauge transformation that induces 
a change of the Lagrangian (59) to 

V = /-ft e «A( e "/2^) 7+7l ^ e «/2 

\ dip 
I dz 

+ (P(e a — - e a i> 1+ ip) + qAe a ip 1+ ip 
\ dip 

Since the first term transforms differently from the other 
terms, the Lagrangian is not time scale invariant. Con- 
sequently, the time parameter must have a physical sig- 
nificance, though that was not imposed in writing the 
Lagrangian. 

We can fix the time scale gauge by requiring p = 
dr/dip = ±1, so proper time becomes our intrinsic time 
variable. That simplifies our model considerably. The 
choice is most conveniently expressed by adopting the 
constraint 



^ = pe iS) = ±1, (89) 

which suppresses the superfluous [3 variable as well. 
There are reasons to associate the minus sign with an 
antiparticle (positron), but we will not explore that pos- 
sibility here and will stick with the plus sign for the rest 
of this paper. 

We could have imposed the constraint (89) at the 
beginning by replacing the spinor ip with the rotor R 
in the Lagrangian. That was not done for two good 
reasons. First, it avoids messing with the constraint 
RR = 1 , which complicates variational derivatives of the 
Lagrangian. Second, the Dirac wave function also has a 
factor pe 1 *, so we want to understand how that relates to 
the similar factor here. Indeed, if the Dirac equation is 
derivable from a superposition of particle paths, as might 
be the case, then the individual proper times must be re- 
placed by a common time variable, so p should emerge 
as a common time scale factor for the ensemble of paths. 

To ascertain the significance of our gauge choice, we 
consider a class of gauge transformations with the form 

V> -> ip' = ipe Fx , (90) 

where, as in (86), the generator T is a constant bivector 
constructed from the vectors 7^. For T = 7271, this is 
the form for an electromagnetic gauge transformation in 
the Dirac equation when it is accompanied by the trans- 
formation of vector potential: 



n n 



(91) 



Inserting this into the Lagrangian (59), we find gauge 
invariance if the following two terms cancel: 



(Vx) • W>'7+V0 = w ' • v x, 



dip' j, 
~chp~ 1+l1 ^ 



^ ^V / 727i7+7iV'' 



dx 



1 I 1 I 



(92) 
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Cancelation occurs if yp'ip'j = 1, which is precisely our 

constraint (89). 

This is an intriguing if not ironic result. Recall that 
Weyl originally introduced gauge transformations as an 
extension of General Relativity to incorporate length 
scale invariance. Subsequently, he renounced that idea 
and introduced the notion of gauge invariance that has 
become standard in quantum mechanics. Here we have 
come full circle to find electromagnetic gauge invariance 
associated with time scale invariance. Perhaps Weyl had 
the right idea in the first place. 



C. Proper equations of motion and first integrals 

Now we are prepared to reformulate our equations of 
motion in their simplest form, with proper time as the 
independent variable. Using the simplification (81), we 
can recast the equation (68) for particle velocity in the 
rescaled form 



u = z = uiei H F ■ u — pu, 



(94) 



where a dynamical zitter frequency has been defined by 
uo = 2m/h. 

The momentum equation (62) can be written 



p = qF ■ u + V$, 



(95) 



where the spin potential for the gradient force is given by 
several equivalent forms: 

$ = $(r,«) = —S(t)-F(z) 

m e 

= —(isu)-F~q\ e F-(uei). (96) 
m e 

For future reference, we note that 

^$ = !^(a r + u .v)$ 

q q 

= (d T S) -F + S-(z- VF) = S ■ F + S • F. (97) 

Finally, from (66) we get the rescaled equation of motion 
for spin 



S = u A p + —F xS-pS, 



(98) 



wherein the coupling with momentum can be simplified 
by (80) to u A p — TO,e 2 eo. 

From the three proper equations of motion (94), (95) 
and (98), we easily derive 



d(pm) d(p • w) 

dr dr 
q d(pS) 



= p • w + p • w 



F + pu- V$ = — (p$), 
m e dr dr 



(99) 



where the identity (42) has been used with (97). Thus 
we obtain an integral of motion for the electron mass: 



The integration constant is readily identified as the free 
electron mass m e , and we recognize $, with its several 
different forms (96), as a variable mass shift due to in- 
teraction. 

To learn more about the relation between mass and 
momentum, we note that the common interaction term 
in cqs. (94) and (95) can be eliminated to get 



dr 



(p — m e pu) = V$ — moj e ei. 



(101) 



This suggests that we define a relative momentum vector 



7r = p — m e pu. 
It then follows that 



n • n — 



d_ 



= TT ■ V$. 



(102) 



(103) 



This describes the rate that spin interaction energy is 
converted to mass. It vanishes when the impressed field 
F is constant, in which case ir 2 is a constant of motion. 

Since m = p-u = mi+m 2 , the various mass parameters 
are related by 

7r 2 = p 2 — 2m e mp = m(mi — m 2 — 2m e p). (104) 

Thus, we have two equations of change (99) and (103) 
for three mass parameters {mi, m 2 , p} or {to, toi— to 2 , p} 
By inspection, no other relations among these parameters 
can be derived from the dynamical equations. This leaves 
us with a gauge freedom to fix one relation among them. 
There are two possibilities with clear physical meaning. 
On the one hand, we could require m\ — m 2 = p 2 — 
(p • u) 2 = to 2 , so there is only one mass parameter. Al- 
ternatively, we can require p = 1 . This appears to be the 
better choice, as it simplifies the equations of motion and, 
as we have already noted, it accommodates gauge invari- 
ance. Here we have learned how it relates to mass pa- 
rameters. However, before committing to a gauge choice, 
there is more to say about the equations of motion. 



D. Rotation dynamics 

So far we have established a system of three equations 
of motion, (94), (95) and (98), for electron velocity, mo- 
mentum and spin. Now it will be shown how these three 
equations can be reduced to a single spinor equation. The 
reduction is important for two reasons: First, it simpli- 
fies solving the equations of motion. Second, it greatly 
facilitates comparison with Dirac theory. 

Reduction is made possible by expressing the momen- 
tum in the form 



p = TOie - to 2 c 2 = RpoR, 
where po = TO470 — to 2 7 2 . 



(105) 



pm = TO e + 



(100) 



Since the velocity is given by u = cq + e 2 = R"f + R and 
the spin is given by S — {%/2)ue\ — (h/2)R-f + -fiR, the 
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evolution of all three observables is determined by the 
rotation dynamics of the comoving frame e M = Rj^R. 
That, in turn, is determined by a single equation of mo- 
tion for the rotor R = R(t) . The normalization condition 
RR = 1 implies that the rotor evolution equation has the 
general form 



it = 1 -sir, 



(106) 



where the bivector rotational velocity fl = O(r) is a spec- 
ified function. It follows immediately that equations of 
motion for the e M are given by 

(107) 



Our problem is therefore to determine the functional 
form of fl. 

Of course, the equations of motion for velocity and spin 
were initially derived from the spinor equation (65). That 
equation is difficult to handle as it stands; however, we 
are now prepared to put it in more tractable form. Ex- 
cluding the duality factor for reasons already explained, 
we can write the spinor in the form ip = p 1 / 2 i?, so (105) 
gives us 

pip = #o- (108) 

Using the identity 707+ = —727+ we have 

Pol+ = (mi + m 2 )7o7+ = "2707+ 

and 7 + po = ^7+70, (109) 

so (65) simplifies to 

fr0'7+ = m-07o7+7i + q\ e Fip^ + . (110) 

The problem now is to solve for tp' . This is complicated 
by the fact that one cannot divide out the null vector. 
The slickest way around this obstacle is as follows. First 
apply (109) to (110) to get 

hip'-f+po = -m 2 7/;7_7i + q\ e Ftpj + p Q . (Ill) 

Then use (108) to get 

= h[(i/jp )' - v^ 7+] = h(p'i> + pi/)' - ipp' h+, 

which, on inserting (110) and the momentum equation 
(95), becomes 

hip'p j+ = Hp(qF • u + V$)V'7+ 

+ PW7+71 + q^eFtp-f+) - VPo7+- (112) 

Adding (110) and (112), we obtain 

2htjj'j + • po = hp(qF • u + V$)V>7+ 

+ ^(P 2 7+ - m 2 7_)7i + q\ e [2p • FV>7+ 

+ FV'bo7+ + 7+Po)] - 7#Po7-H 
This can be solved easily for ip' using 



and further simplified with 

p 2 7 + - to 2 7_ = m[(mi - m 2 )"f+ - (mi + m 2 )7_] 
= mpo7o72, 
Pol+ = ( m i1o - m 2 7 2 )7+ = m'7o7+. (113) 

Finally, we obtain the desired spinor equation: 



v 2 



-pe e 2 ei 
n 



^F+ 



— ( /'V«I> — F • n — m c,i ] 11 

m \ 



m e 



where it = p — m e w = p — m e pu, as before. 
Expressing (114) in the general form, 

V'V -1 = P '/p + n = p + n, 



V>, (H4) 



(115) 



we can separately equate scalar and bivector parts. The 
scalar part gives us 



-F-{pAu), (116) 



p' m' (pm)' p „ ^ 
- — I = - — - = —u • V$ - - 

p m pm m 1 

which is identical to the equation (98) that gave us the 
mass integral of motion (99). The bivector part gives us 

2 q 

O = ^rpe e 2 ei H F 

n m e 

1 / a \ m' 

H —F * 7r I Au -\ e 2 e . (117) 

m \ ra e J m 

This is the general expression for rotational velocity that 
we were seeking. It holds for any mass gauge choice, 
including p = 1, which we now adopt for the balance of 
this paper. 

The first term in (117) can be interpreted as a kinetic 
rotational velocity: 

2 2 

^fe = ^pe e 2 ei = ^( m i e 2 ~ m 2e a )ei- (118) 

It is a spacelike vector (since ill = —Ap 2 /h 2 ), so it gener- 
ates a spacelike internal zitter rotation, as is evident from 
its contribution to the velocity equation: • u = uiei, 
where the zitter frequency lu = 2m jh, is determined by 
the mass equation m = m e + $. The second term 
(q/m e )F can be regarded as a generalization of the 
Lorentz force, as it determines extrinsic bending of the 
particle history due to external fields (in (94)), along with 
a torque on the spin (in (98)). The remaining terms ac- 
count for shift in the internal rotation rate to accommo- 
date change in mass, but they do not directly contribute 
to the equations for velocity and spin. 

Further insight into the rotational velocity is obtained 
by multiplying (114) by the spin bivector (76). The m 
term cancels out (when p = 0), so only the first two terms 
survive to give us the elegant relation 



7+ - Po = 2(7+Po+Po7+) =p-u = m, 



ns = n- s + nx s + nAS = - P u + -^-fs. (119) 
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Separation into parts of homogenous grade yields: 

q 



n-S = -m+-^-F-S 



(120) 



nxS = uAp+—F x S, 



Q, A S = —F A S. 



(121) 



(122) 



Here we see that the mass integral of motion m — $ = 
m e = — fl-S can be interpreted as constraining projection 
of the spin onto the rotational velocity to be a constant 
of motion. 

It is noteworthy that (with p = 1) the mass integral 
of motion can be reformulated as an expression for the 
zitter frequency: 



+ -L F . = ,123) 

m„ dr 



where, anticipating further analysis, a zitter phase an- 
gle Lp has been introduced. We shall see that this an- 
gle is completely analogous to the phase angle in the 
Dirac equation and, a fortiori, to the phase angle in 
Schroedinger's equation. Note that the phase angle may 
depend on the electromagnetic vector potential, but the 
frequency shift is gauge invariant, as it has the form of 
a dynamical flux integral on the circulating particle his- 
tory. The present context also suggests that the phase 
ip may be a more fundamental time variable than the 
proper time. 

Generalizing the free particle case (58), we identify 



A 



dr_ 

dip ' 



(124) 



as the variable zitter radius. This interpretation is sup- 
ported by the form of equation (94), wherein A can be 
identified as an intrinsic radius of curvature for a helical 
particle history. 

The first curvature k x of a particle history measures 
the bending rate orthogonal to its velocity. In the particle 
equation of motion (94) the direction of helical bending 
is given by the zitter vector e\. Hence, the first curvature 
is given by 



-u • e\ 



-?-F.(uei) 



A- 1 . (125) 



According to (101) this quantity is rigorously constant! 
However, it is composed of two parts: the first can be 
regarded as an intrinsic curvature while the second is an 
"extrinsic curvature" due to external forces. This tells 
us that in response to external forces the intrinsic cur- 
vature is adjusted to maintain an overall constant value. 
Its inverse is the curve's radius of curvature, which we 
recognize as the free particle zitter radius. 



VII. ANALYSIS AND IMPLICATIONS 

The formulation of our model for a lightlike electron 
with zitter is now complete. In principle, the model can 
be applied to any problem in quantum mechanics. The 
present discussion is limited to a general analysis of so- 
lution techniques and a survey of promising possibilities 
for applications and experimental test of the model. 



A. Zitter center 

An important unresolved issue is how best to define the 
relation between the timclikc reference history x = x(t) 
and the lightlike particle history z — z(t). Relative to 
the instantaneous rest frame defined by v = x, at any 
time r the two histories are separated by a displacement 
vector 



x. 



(126) 



If this is defined so that \r\ ~ A e , we can regard the refer- 
ence history as specifying a center of curvature, around 
which the helical particle revolves as time progresses. An 
essential condition on a zitter center history x{t) is 



x • u 



x • z 



(127) 



as that defines the proper time. What we need is a 
tractable equation of motion for x that satisfies these 
constraints. This is akin to defining an equation of mo- 
tion for the guiding center of a helical orbit in plasma 
physics. Two possibilities present themselves immedi- 
ately: x = p/nip or x — e . Unfortunately, their equa- 
tions of motion imply that they both "wobble" with the 
zitter frequency, and if we try to smooth that out (be- 
low), there is no guarantee that the history x(t) cannot 
drift outside the helix z(t). 

An alternative with attractive physical interpretation 
is to specify x(t) by defining a zitter radius vector r = 
— Aei with variable length specified by (124). However, 
since the first curvature k\ is a constant of the motion, 
a better choice may be r e = — A e ei, defining a zitter 
radius vector of fixed length. The felicity of that choice 
is confirmed by noting that it gives the spin potential 
(96) the perspicuous form 



$ = qF • (u A r e ). 



(128) 



As there is no crucial physical issue involved, we can de- 
fer choice among the alternatives to be decided by conve- 
nience in specific applications. The best choice may well 
depend on what approximations we adopt. 



B. Zitter averages and approximations 

As zitter fluctuations are so rapid, it is most conve- 
nient to separate them from zitter means, which are more 
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directly observable. The velocity v — x defines an instan- 
taneous rest frame for the electron at each time r, so we 
define the zitter mean as an average over the free particle 
zitter period that keeps the zitter center velocity v and 
the spin vector s = (h/2)es fixed. With an overline to 
denote average value, basic zitter means are specified by 
the expressions 

v = v = x = u = e , ei=0 = e 2 , e 3 = e 3 . (129) 

Consequently, the mean of the spin bivector is 



S = —uei = — e 2 ei = isv . 
2 2 



(130) 



This approximation ignores variations in zitter radius 
and mass over a zitter period. 

Since the electromagnetic field F(z) — F(x + r) acts 
at the location of the particle, to get its effect on the 
zitter center, we expand with respect to the zitter radius 
vector. Accordingly, the average field at the zitter center 
is given by 



F(z) = F{x) + r • VF(x) + - (r • V) 2 F(x) + . . . . (131) 

Of course, r — 0, so to first order F = F{x). 

Neglecting second order effects gives us the mean mass 

m = m e + $, with $ = —S-F(x). (132) 

m e 

The mass shift $ should be recognized as a generalization 
of the Zccman interaction in atomic physics. 

In the present approximation the mean momentum 
aligns with the zitter center velocity: 



p = TOie - m 2 e 2 = m e v, 
where we have used 



mi = m e , so to 2 



m 



mi = <I>. 



(133) 



(134) 



Consequently, the momentum equation (95) can serve as 
an equation of motion for the zitter center. Projecting 
out the effect of the mass derivative and averaging, we 
obtain from (95) a mean equation 



m e v = qF • v + v • (v A V)3>. 



(135) 



This will be recognized as the classical equation for a 
charged particle with intrinsic spin, so it can be regarded 
as the classical limit of the zitter model. With F = F(x), 
the first term on the right is the classical Lorentz force, 
while the second is the Stern-Gerlach force. 

For slowly varying electromagnetic fields there is not 
much difference between equations (95) and (135), save 
that replacement of the Lorentz force qF • v by qF • u 
produces a wobble of v about its mean value. The big 
difference in the equations comes with electromagnetic 
fields oscillating with a frequency close to the zitter fre- 
quency, for then resonance can occur, as we shall see. 



Our average model is not complete without an equation 
of motion for the spin. We can get that by computing 
the mean rotational velocity from (116), with the result 



Q = eie 2 w e - — v A V$ + —F 
m e m P 



/ 2$ $ \ 

eo ei— e 2 — . (136) 

\ n m 



The last pair of terms might be regarded as a second 
order effect, but if the zitter vectors ei,e 2 are "dotted 
out" it surely becomes first order. The mean motion of 
the comoving frame is now given by 



— ^ * e^ , 



(137) 



which determines the mean spin motion as well. Let's 
call this the minimal zitter model, because it is the sim- 
plest model that captures the main features of zitterbe- 
wegung. Obviously, the model can be extended by in- 
cluding higher order moments in the averaging. 

In the next Section we see that the minimal model 
corresponds closely to Dirac theory, with v = u corre- 
sponding to the direction of the Dirac current, while S 
matches the Dirac spin. That may be surprising, since 
the spin bivector here satisfies the null constraint S 2 — 0, 
and no such thing appears in standard accounts of Dirac 
theory. The discrepancy is resolved by considering the 
spacetime split [1]: 



S = 



h h 

-ue\ — —isu 



2 2 
= (S • v + S A v)v 



Sv 2 

= (m e r + is)v . 



(138) 



where the zitter radius vector and mean spin are given 

by 



and 



m e r = ei = S • v 

2 



S = (S A v)v = isv , 



(139) 



(140) 



with the sign chosen so the radius vector points away 
from the zitter center. 

This is a split of the spin S into a rapidly rotating 
part m e rv and a slowly precessing part S. The specifi- 
cation of spin vector here is algebraically identical to the 
one in Dirac theory [1] . The rapidly rotating unit vector 
f = — ei is also inherent in Dirac theory, but its physical 
significance has been overlooked. In Dirac theory, identi- 
fication of the radius vector and the spin vector as parts 
of the null bivector S has not been made heretofore. The 
implication is that the spin observable in Dirac theory is 
only the zitter average S of the full spin S. Nevertheless, 
zitter is still present in the Dirac equation as a rotating 
phase factor in the wave function, and we shall see that 
there are indeed possibilities to detect it. 
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The "v-split" of the spin in (138) entails a w-split of 
the electromagnetic dipole moment into electric and mag- 
netic parts: 



— S = d v + ifi v 



with 



d v 



-qrv, 



ii v = — sv 



(141) 



Note that the electric dipole has constant magnitude 
±|d„| = qX e = qh/2m e , so it only rotates. A correspond- 
ing n-split of the external electromagnetic field has the 
form F = E„ + iB v , where electric and magnetic parts 
are given by 



E v = -(F-vFv), 



iB v = -(F + vFv). (142) 



Of course, these are electric and magnetic fields as "seen" 
in the instantaneous rest frame of the zitter center, not 
to be confused with fields in a "lab frame." 

The induced mass shift can now be expressed in the 
physically perspicuous form 

q 



-S • F — d v • E„ — ij, v • B v . 



(143) 



This is a general result, holding exactly in every applica- 
tion of the zitter model. As d v is a rotating vector, this 
result is consistent with the well established experimen- 
tal fact that the electron does not have a static electric 
dipole moment. Of course, if the E„ field is slowly vary- 
ing, the effective shift is reduced to the Zeeman interac- 
tion $ = — /j, v • B v . 

However, the spin split described here has a draw- 
back, namely, that the reference direction v is constantly 
changing with time, so it is difficult to compare spin di- 
rections at different times. We show next how that draw- 
back can be eliminated by examining the spin split more 
thoroughly from a different point of view. 



C. Zitter and spin in the electron rest frame 

Transformation of the time dependent instantaneous 
rest frame specified by velocity v = v(t) to an inertial 
rest frame specified by a constant vector 70 is a boost, 
specified by 



v = L"f a L = L 2 "f or 
L 2 = vjo = v Q (l + v), 



(144) 



where vo — v • 70 is the relativistic time dilation factor 
relating particle proper time r to inertial time t, that is, 



1 



v = v • 70 = 



(1 



dt 



~dt 



(145) 



is the relative velocity of the particle in the inertial refer- 
ence frame. Let us refer to this frame as the electron rest 
frame, though we shall generalize the definition some- 
what in a later subsection 



We can solve (144) to get an explicit expression for the 
rotor L that generates the boost [1] : 

L=(vy )2= r = r . (146) 

[2(1 + V )]2 [2(1 + V )}2 

However, it is usually easier to work with the simpler 
form for L 2 in (144). 

The deboost of the u-split for spin S in (138) to a split 
in the electron rest frame is now given by 

S = LSqL where So = —m e r + is, (147) 
with electric and magnetic moments defined by 

— So = -qr + i— s = d + ifi. (148) 

m e m e 

The deboost of the particle velocity has the form 

uv = L(l + u)L. (149) 

Since spin and velocity are related by Su = 0, from (147) 
and (149) we obtain 

5 (1 + u) = (-m e r + is)(l + u) = 0. (150) 

Separating parts of homogeneous grade, we ascertain 

r-u = 0, s-u = (151) 

and 

s = —im e YVL = m e r X u. (152) 

This makes spin look like a classical orbital angular mo- 
mentum. However, the velocity u = z is not generally 
the derivative of the radius vector r, as we shall see. 

It is often useful to express the relation of spin to zitter 
in terms of a rotating orthonormal frame: 



e fc = Le k vL = Uer k U, 



(153) 



where ei = r, e 2 — u and e 3 = s. However, it is phys- 
ically more enlightening to first study the dynamics of 
spin and zitter more directly. 

To treat interaction with external fields in the electron 
rest frame, we need the deboost of an electromagnetic 
field: 



LFL = E 







iB , 



where 



F = E + iB = LF L 



(154) 



(155) 



expresses the field in terms of electric and magnetic fields 
defined in the lab frame. 
Note that 

< SF >=< LSqLF >=< S LFL > 

= < S F >=< (d + i/x)(E + iB ) > . (156) 

Whence the spin potential (= mass shift) assumes the 
form 

$ - —S ■ F = -2-So • F = d ■ E - • B . (157) 

m e m e 
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However, this is not the final form required for physical 
interpretation. 

We still need explicit expressions for the deboosted 
fields Eo and Bo in terms of lab fields E and B. Us- 
ing (144) we have, with obvious notation, 



F = LFL = F|| +F±L 2 

where 

F|| v = vF|| , Fj_v = -vFi. 
Whence 

E = E|| + D Ei + v v X B, 

where 

E|| = (E-v)v, E_l = E-E|| 
Similarly, 

B = Bii + v Q B ± - v Q v X E. 



(158) 
(159) 
(160) 
(161) 
(162) 



Insertion into (157) gives explicit expressions for the in- 
teraction of the electron's electric and magnetic dipoles 
with external fields. 

Physical interpretation of spin dynamics is facilitated 
by transforming the spin equation of motion (98) to the 
instantaneous rest frame. Thus, 

S = L(S + n Q x So) = L(me 2 + —F x S )L, (163) 

where the rotational velocity £Iq is determined by differ- 
entiating the boost (146), with the result [1]: 



LQ L = Q v = 2LL = 



v A (v + 70) 

T+V~o 



v v .viv X V 

«ov + — + • 

1 + vo l + vo 



(164) 



As is evident in equation (163), this quantity acts like an 
effective electromagnetic field induced by acceleration of 
the rest frame. The magneticlike term at the right side 
of the equation is responsible for the classical Thomas 
precession of the spin. 

Inserting the split (147) with (148) into (163), we have 



So = m e r + is 



m e 



= me 2 + {F -— -fi ) x (d + i/i). (165) 
Now we introduce the split 

Tfl — - 711 

& + ib = F Q --j-n = L(F--j-n v )L, (166) 

noting that expressions for a and b can be obtained from 
(164) by the deboost (158). Using 

(a + ih) x (d + in) 

= aAd-bA/i + i(aA/i + bAd) (167) 



with a A d = i(a X d), we split (165) into coupled equa- 
tions of motion for zitter and spin: 



m e r = 77ie 2 + /xXa + dxb, 
s = axd + fixb. 



(168) 
(169) 



These equations are helpful for analyzing the dynamical 
behavior of spin and zitter vectors. 

The spin equation (169) is most familiar. Its last term 
H X b is recognized as the usual spin precession torque, 
including the Thomas precession, as already noted. The 
term a X d is something new. First note that 

s • s = (a X d) • s = (d X s) • a = = e 2 • a. (170) 

To prove that this quantity does indeed vanish, we use 
(166): 

ma • e 2 =< LaLu Ap >~< (F -Cl v )u Ap > 

[u-v + v-u]=0. (171) 



It follows from this constraint that we can write 
a X d = s X d (s • a). 



(172) 



This torque has several interesting properties. First, it 
rotates s about the vector e 2 , which is the direction of 
zitter circulation in the rest frame. Second, the torque 
decreases in magnitude until it vanishes at s • a = 0. 
Third, its zitter average is zero if a is slowly varying. 
These properties suggest that the term is a good candi- 
date for a spin polarization torque. 

The zitter equation (168) has similar properties, of 
course, with the additional term me 2 expressing the high 
frequency zitter rotation. We only note that the fi X a 
torque ensures that r rotates around e 2 along with s. 



D. Zitter model for a static potential 

The interaction of atoms and crystals with an electron 
is usually modeled with a static potential. For that pur- 
pose we introduce a static vector potential defined in the 
incrtial frame of 70 at each time t and position x by 

qA = Vj where V = V(x). (173) 

This determines an electromagnetic field 

qF = VU A 70 = -70 VU = -70 A W = - W. (174) 

This, in turn, determines an electromagnetic force on the 
particle 

qF -u = Wu • 70 - 70U, where V = u • W. (175) 

Inserting this into the momentum equation (95) we ob- 
tain 



d_ 



(P + U70) = -V(Uu- 7 o + $). 



(176) 
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Next, we introduce the space-time split [1] pr/o = p • 70 + 
p A 70 = Po + p and separate spatial and temporal parts. 
This yields a conserved energy: 

E = Po + V (177) 

and, after factoring out v n w u • 70, a momentum force 
law 

^ = -V(T/ + V 1 $), (178) 

where V70 = «o(l + v ) with v = dx/dt and p = mvov. 

We still need an explicit form for the spin potential. 
From (175) we get 

$ = X e qF • (uei) = m • 70 ei • W - 70 • e 1 V. (179) 

However, it is better to use (160) and (162) in (157) to 
get the form 

$ = q\ e [e • (E|| + u E ± ) + s • (w v X E)] 

with qE = -W, (180) 

where e denotes the unit zitter vector. This has the ad- 
vantage of clearly separating the high frequency zitter 
from the low frequency spin contributions, and it ex- 
hibits the zitter radius A e as governing the strength of 
the interaction. 

These results can be applied to experimental search for 
observable effects of electron zitter. As the energy conser- 
vation law (177) is identical to the usual one, zitter will 
be manifested only in momentum fluctuations, which are 
so small and rapid that they are observable only in res- 
onances. Consider Mott scattering by a Coulomb field, 
for example. In low and high energy ranges the cross 
section will not be significantly affected by the very high 
frequency zitter rotation, so the standard result should 
be obtained [23]. However, in an intermediate range 
where the impact parameter is on the order of a zitter 
wavelength, the zitter structure of the particle trajectory 
should have a significant effect on the scattering. How 
big an effect awaits calculation! 

As explained below, most resonances involve zitter field 
theory, which is beyond the scope of the present treat- 
ment. However, a new possibility for amplifying zitter 
resonances in crystals has recently appeared in electron 
channeling experiments, to which we now turn. 

E. Zitter resonance in electron channeling 

When a beam of electrons is channeled along a crystal 
axis, each electron is subject to periodic impulses from 
atoms along the axis. When the energy of the beam is 
adjusted so the crystal period matches the period of the 
electron's zitter dipole, a resonant interaction may be ex- 
pected to alter the distribution of transmitted electrons. 
Indeed, a pioneering experiment in search of such a res- 
onance has already been performed, but without knowl- 
edge of the dipole interaction mechanism described here. 



The anticipated resonance was observed at close to the 
de Broglie frequency, which is precisely half the zitter 
frequency [24] . Our purpose here is to show how this re- 
sult can be explained quantitatively by the zitter model. 
This confluence of theory and experiment provides ample 
grounds for repeating the experiment with greater accu- 
racy to confirm the results and look for details suggested 
by the theory. 



1. Experimental specifications 

The anticipated resonant energy is easily calculated 
from the de Broglie's (circular) frequency ujb = m e c 2 /h. 
One of de Broglie's original insights was that the fre- 
quency of a moving electron observed in a laboratory 
will be u>l = wb/7, where 7 = vq is the relativistic time 
dilation factor. The distance traversed during a clock 
period is d = 2wcf3/ujL = hp/(m e c) 2 . For the silicon 
crystal used in the experiment, the interatomic distance 
along the < 110 > direction is d — 3.84 A, which implies 
a resonant momentum p = 80.874 MeV/c. 

In channeling the maximum crystal potential is a few 
hundred electron volts at most, so in the 80 MeV region 
of interest to us, the effective electron mass M = jm e = 
E/c 2 is constant to an accuracy of 10 -5 , and 7 = 158. 

In axial channeling electrons are trapped in orbits spi- 
raling around a crystal axis. To a first approximation, 
the crystal potential can be modeled as the potential for 
a chain of atoms, so it has the form 

V(x) = V(r, z) = U(r)P(2nz/d), (181) 

where x(i) = r + zer z is the particle position from the 
first atom in the chain, with r = |r|. The longitudinal 
potential P(2irz/d) = P(uot) is periodic with a tunable 
frequency ujq = 2irz/d that varies with the energy E. 
Note that at the expected resonance z = d/T[ Jl so uj — 
2it/Tl = ljbH corresponds to the de Broglie frequency. 

Our problem is to calculate perturbations on the trans- 
verse component of the momentum vector, as that can 
remove electrons from stable orbits in the beam. The 
transverse component of equation (178) has the familiar 
form of a nonrelativistic equation: 

Mr = -r(Pd r U + 7~ 1 <9 r $), (182) 

while the longitudinal component has the form 

Mz = -{Ud z P + 7 _1 5 z $), (183) 

where now the overdot indicates differentiation with re- 
spect to "lab time" t. In the energy range of interest, it 
is easy to show that oscillations in the transverse veloc- 
ity can be ignored, so we regard i as a constant tunable 
velocity close to the speed of light as already assumed 
above. 
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2. Crystal potential 

To proceed further, we need an explicit model of the 
crystal potential. For analytic simplicity, we approximate 
the potential by the first two terms in a Fourier expansion 
with respect to the reciprocal lattice vector. So write 
(181) in the more specific form 



V = U(r)(l + cosw i), 



(184) 



where the coefficient of the second term is set to make 
the potential vanish between atoms. The first term is 
the potential for a uniformly charged string, which (in 
its simplest version) has Lindhard's form [25, 26]: 



U(r) = -fcln[l + (Ca/r) 1 



(185) 



where Ca — 0.329A, and (with Z = 14 for silicon) 
the coupling constant k — Ze 2 /d — 52.5 eV is the 
product of electron charge with charge per unit length 
of the string. The constant a = 0.885Z~ 1//3 ao = 
(0.885)(14)- 1 / 3 (0.529A) = .190A is the Fermi-Thomas 
screening radius, and the constant C 2 = 3 is a fairly 
accurate fit over the range of interest. 

The string potential is defined as an average over 
atomic potentials; thus, 



U(r) 



-f 

J — < 



Kt om [(r 2 + z 2 ) 1/2 ] 



dz 



(186) 



where the screened atomic potential is given by 

Ktom(fl) = -^-<p(R/a). (187) 
To get (185), Lindhard used the screening function 

<p(R/a) = 1- [1 + iCa/R) 2 ]- 1 / 2 . (188) 

Actually, the experiment is not very sensitive to the 
shape of the potential, so long as it is sufficient to bind 
the electron to an orbit around the crystal axis. In the 
first approximation V = U(r) the projection of the elec- 
tron orbit onto a transverse plane looks like a precessing 
ellipse or rosette. The second term U cos Lj t is ignored 
in most channeling calculations, as it merely produces 
small harmonic oscillations on the radius of the precess- 
ing ellipse. However, the periodicity of the second term 
is essential for resonance in the zitter perturbation of in- 
terest here. 



3. Classical channeling orbits 

Ignoring the zitter perturbation in (182) for the time 
being, we seek to ascertain the effect of the periodic factor 
P(uiot) = 1 + cosu i on the orbital motion. With an 
obvious change of notation, we can represent the radius 
vector in the complex form r = re lS , where the imaginary 



i is the bivector generator of rotations in the transverse 
plane. Then equation (182) assumes the complex form 

TJI 

Mr = -U'Pe ie = - — (1 + cosw t)r, (189) 

with U' = d r U. 

We are interested only in radial oscillations, so we use 
conservation of angular momentum L = Mr 2 9 to sepa- 
rate out the rotational motion. With the periodic driving 
factor omitted, equation (189) admits the energy integral 

E± = ^Mr 2 + W(r) 



where 



W(r) 



2Mr 2 



U(r). 



(190) 



Let us expand this around a circular orbit of radius r , 
and for quantitative estimates take r = 0.50A as a rep- 
resentative intermediate radius. For r = re lS °, equation 
(189) gives us 



9o 2 = 



3. 

Mr n 



L 2 



M 2 r\' 



where C/g 



U'(ro). In terms of x 
(190) gives us 



Ej 



1 



Mx 2 + W 



:W'x 2 , 



(191) 

ro, expansion of 
(192) 



where 



W = U +r U'j2, W^ = U^-U 
Wg = US + 3^/ro, 



= 0, 



with U = -18.9 eV, 
9k 

u> = drU = fl 



(Ca/r) 



l + (Ca/r) 2 
with r U^ = 31.7 eV, and 
"3 + (Ca/r) 2 



u» = - u - 

r 



1 + {Ca/rf 



(193) 



(194) 



(195) 



with r 2 U{l = -76.0 eV. 

Differentiating (192) and reinserting the periodic driv- 
ing factor, we obtain the desired equation for radial os- 
cillations: 



x + fl (l + COSLOot)x = 0, 
where, for mass M at the expected resonance, 

roUg 



(196) 



n 2 ^ vr ° 



M 



3U 



Mr 



3 x 31.7-76.0 c 



3 ). ( 197 ) 



80.9 x 10 6 
' 1 . This should be compared with 



so O q = 4.21 x 10 15 s 
the O = 4.75 x 10 15 s _1 from (191) and the expected 
resonant frequency ujq = uob/i = 4.91 x 10 18 s _1 . We 
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note that the distance traveled in one orbital revolution 
is d r = 2wc/e = 3.97 x 10 3 A = 0.397^m. Thus, the 
orbit makes 2.52 revolutions in passing through the one 
micron crystal, so the orbital revolutions are only weakly 
coupled to the high frequency radial oscillations. 

Equation (196) is a special case of Mathieu's equation: 

x + q(l + hcosLut)x = 0. (198) 

According to Floquet's Theorem, this equation has solu- 
tions of the general form [27, 28] 

oo 

x{t) = e st a n e lnuJt . (199) 

n— — oo 

Substitution into (198) yields a three term recursion rela- 
tion that can be solved for s and ratios of the Fourier co- 
efficients a n . In general, the Floquet exponent s = s(q, h) 
is a complex constant, so it determines whether solutions 
are unstable or periodic. 

Since £1 << uj in equation (196), its Floquet expo- 
nent s = ifl is pure imaginary, and its recursion relations 
give 

3 

J^ = J^>_ = C% forn>l. (200) 
a n -\ a-n+i 2n 2 Lu£ 

Therefore, the solution is dominated by the first order 
term, with the particular form 

x(t) = a(cosQt) cosuj t = ^(cosw + t + cosw_i), 

w± = ujq ± n. (201) 

As might have been anticipated, this describes a har- 
monic oscillator with high frequency loq and a slowly 
varying amplitude with frequency fi, which is equivalent 
to a sum of two oscillators with frequencies u>± separated 
by uj + — uj- = 2fl. 

We shall see that, at the resonant frequency uo = 
Wb/7, the frequency shift = (0.857 x 10~ 3 )oj is the 
right order of magnitude to contribute to experimental 
effects. Moreover, this quantity has been estimated at 
the particular radius r = 0.50A, and it may be larger by 
an order of magnitude for smaller radii of experimental 
relevance. Accordingly, a distribution of il values will 
contribute to the experiment. 

A limit on the maximum radius of a channeled elec- 
tron is set by the requirement that W(r) in (190) must 
be negative for a bound orbit. A sharper limit is set by 
neighboring atoms. The total crystal potential can be 
modeled as a sum of chain potentials isomorphic to the 
one for the channeling axis. For larger radii perturbation 
from other chains can induce transition to a neighbor- 
ing chain, with the result that channeled electrons jump 
randomly from chain to chain. We are interested in this 
effect only to the extent that it affects the distribution 
of electrons transmitted by a single channel. However, 



a quantitative estimate of such transitions will not be 
attempted here. 

In silicon, the closest chain to a < 110 > channel is 
at a distance of 1.36 A with atoms located at positions 
z n+i/2 = (n+ l/2)d alternating with the positions z n — 
nd along the channel. Considering the slow precession of 
a channeled orbit at the resonant frequency, this chain 
will resonate with it for hundreds of atomic steps. In 
fact, the interaction might lock onto the orbit to prevent 
precession during resonance. This may indeed have a 
substantial effect on channeled electrons, but we will not 
investigate it further here. 

4- Zitter perturbations 

Now we are prepared to consider the effect of zitter 
perturbations on the orbit. To evaluate the zitter poten- 
tial we use equation (180) with E|| = (E • v)v and 
v « cr z = Vz, and we ignore the spin term (because it 
vanishes on averaging over spin directions). By the way, 
we don't need to use the spin equation of motion (98) 
in our calculations; we only need the fact that it implies 
that the unit zitter vector e rotates rapidly in a plane 
that precesses slowly with the spin s. 

Inserting the electric potential into (180), we get 

$ = -\ e [je-rU'P + e-a z UP'}. (202) 

Since r is effectively constant over a zitter period, we can 
make the zitter oscillations explicit by writing e • f = 
cos(w ztf '7 + 5) and e • er z = sm(ujzt/j + 5'), where uiz is 
the zitter frequency. Thus, 

$ = -A e [7f7'Pcos(w z i/7 + s ) 

+ UP'sin(wzth + 5% (203) 

where the smaller second term can be neglected. 

The spin potential has a twofold effect on electron mo- 
tion: first, as a shift in zitter frequency according to equa- 
tion (123); second, as a perturbation of the momentum 
in equation (182). For the frequency shift we get the 
explicit expression 

-fU' 

ujz = oj e cos(w e i/7 + S). (204) 

m e c 

The modulus of the oscillating term has the esti- 
mated mean value ~/U'c/m e c 2 = (158)(31.7eV/0.5A)(3 x 
10 18 As-V(0.511 x 10 6 eV) = 1.96 x lO^s" 1 . Compared 
to uj e = 1.55 x 10 21 s _1 , this quantity is too small by 10~ 3 
to play a role in the present experiment. However, the 
value of V, estimated here at a radius of r = 0.5A, 
may be regarded as lower bound, as the logarithmic po- 
tential is a very soft approximation to a realistic po- 
tential, which is much sharper at the screening radius 
a = 0.190A where the logarithmic potential is invalid. 
We shall see that close resonant encounters with atomic 
nuclei play a dominant role in the experiment, so the os- 
cillating term might indeed contribute a mean frequency 



18 



shift on the order of 10 _2 o; e , which could show up in ex- 
perimental data. Replacement of the soft periodic mod- 
ulus P(u)ot) = (1 + cosw t) by one that is more sharply 
peaked at the atomic sites may strengthen this conclu- 
sion. That being said, for the rest of our analysis, it 
suffices to assume iqz — w e . 

Inserting (203) into (182), we get an explicit expression 
for the zitter perturbation term: 

-7 _1 (V$)_l = -f-Hdr® 

= X e rU"Pcos(uJzt/"/ + 5). (205) 

5. Zitter resonance 

The task remains to show that this interaction can 
produce a resonant amplification of the electron's orbit. 
Inserting it into equation (182) with a convenient choice 
of phase and writing to — coz/j, we get the equation of 
motion 

d 2 r U' A 

M — = -r— (1 + cos w t) (1 + coswf), (206) 
at z r R 

where R = -U'/U" = r[l + (Ca/r) 2 }/[3 + (Ca/r) 2 } is 
an effective (screened) radius. Ignoring the amplitude 
modulation as determined in (201), we can reduce this 
to a radial equation 

x + u%(l + ^costot)x = 0. (207) 
R 

Of course, we can replace coo in this equation by ui± to 
get two separate resonant peaks. 

For h = X e /R constant, (207) is Mathieu's equation 

(198) , so let us solve it using h = X e /Ro = 1.931 x 
10~ 3 /0.208 = 9.283 x 10~ 3 as an approximation. For 
small values of h such as this, equation (207) has a first 
order resonance at to = 2uj . An easy way to see that is to 
regard (207) as an equation for a driven harmonic oscil- 
lator with driving force proportional to its amplitude. In 
the experiment, the "natural frequency" loq was adjusted 
by continuously varying the initial energy of the electron 
until a resonance was found. As for a periodically driven 
harmonic oscillator, resonance occurs when the driving 
frequency equals the natural frequency. In this case, that 
means to — loq =w orojo = to/ 2 = ujz/2j, which explains 
why the resonance occurs at half the zitter frequency — 
a surprising result until contemplation shows that it is 
an obvious feature of parametric resonance! 

More generally, it can be shown that (207) has reso- 
nances at 2uj = nto, for n = 1,2,.... We demonstrate 
that explicitly for the first order resonance by truncating 
the series in (199) to the form [29] 

rfr tot , . Lot. ,„„„s 

x = e s *[acosy +6siny]. (208) 

For a resonant solution, the Floquet exponent s must 
be real and positive. [Actually, to get this form from 

(199) we should use the Floquet coefficient s — iuo/2 and 



incorporate its imaginary part in the series. Thus, there 
is some ambiguity in the definition of Floquet exponent.] 
To validate the trial solution and evaluate its param- 
eters, we insert it into the differential equation. Using 
trigonometric identities such as 2 cos A cos 2 A = cos A + 
cos3A to reduce products to sums and dropping higher 
order terms, we obtain 

stab + [s s +lu 2 - (uj/2) 2 + (h/2)u 2 ]a = 0, 
sua - [s s + (lo/2) 2 - (h/2)tol]b = 0. 

For a near resonance solution we write to = 2uj + e and 
neglect second order terms to get 

sb-^[e- ^huj o ]a = 0, 

sa+ i[e+ ^hu) ]b = 0. 

These equations can be solved for the coefficients pro- 
vided 

s 2 = 1 -[(lhco ) 2 -e 2 ]> 0. (209) 

Thus, we have a resonance with width 

Alo = 2e = huj = 9.283 x 10~ 3 w , (210) 

or 

Ap = /i(80.9 MeV/c) = 0.751 MeV/c. (211) 
And for the amplification factor at resonance we have 

1 7T 

st = -hcoot = -hn = 1.46 x 10~ 2 n, (212) 

where n is the number of atoms traversed in a resonant 
state. Since In 2 = 0.693, this implies that the amplitude 
is doubled in traversing about 50 atoms. 

Now note that the value of h = A e /i?o used in (210) 
and (212) applies to only the subclass of orbits for which 
ro = 0.50A. For smaller radii the values can be much 
larger. In principle, to get the width of the ensemble 
of orbits we should replace h in (210) by its average h 
over the ensemble. However, the result will probably not 
differ much from the typical value we have chosen. 

Similarly, the amplitude factor in (212) will have a 
distribution of values, and the doubling factor will be 
reached much faster for orbits with smaller radii. Pre- 
sumably, random perturbations (such as thermal fluctua- 
tions of the nuclei) will limit the resonant state coherence 
length to some mean value n. Consequently, states with 
smaller ro will be preferentially ejected from the beam. 

Even more to the point, the perturbation parameter 
h = \ e /R is not constant as we tentatively supposed but 
increases rapidly as the electron approaches a nucleus. 
In resonance the value of h close to each nucleus domi- 
nates the perturbation, so its effective mean value is much 
smaller than the estimate for constant r . Evidently, res- 
onant interaction may eject electrons with small ro in just 
a few atomic steps. 
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6. Experimental implications 

The predicted resonance width in (211) is in fair agree- 
ment with the width in the channeling experiment data 
[24] , considering uncertainties in the value of h and such 
factors as thermal vibrations that may contribute to 
damping. Damping can only narrow the width and de- 
stroy the resonance if it is too severe. 

We need to explain how the orbital resonance is man- 
ifested in the experimental measurements. Two scintil- 
lators, SC2 and SC3, were employed to detect the trans- 
mitted electrons. The larger detector SC3, with a radius 
about 3 times that of SC2, served as a monitor while the 
smaller detector served as a counter for a central por- 
tion of the beam. The measured quantity was the ratio 
of SC2 to SC3 counts as momentum was varied in small 
steps (0.083%) over a 2% range centered at the expected 
resonance momentum 80.9 McV. An 8% dip was observed 
at 81.1 MeV. Orbital resonances may contribute to this 
effect in at least two ways: first, and perhaps most impor- 
tant, by increasing the probability of close encounter with 
a nucleus that will scatter the electron out of the beam; 
second, by increasing the duration of eccentric orbits out- 
side the central region. Overall, resonant interactions will 
be strongest on electrons confined to the central region. 
These observations suffice for a qualitative explanation 
of the observed dip. A quantitative calculation will not 
be attempted here. 

The most problematic feature of the experiment is the 
0.226 MeV/c difference between observed and predicted 
resonance energies. If estimation of the experimental er- 
ror was overly pessimistic, that indicates a physical fre- 
quency shift. The most likely origin for such a shift is the 
frequency split in (201). The experiment was not suffi- 
ciently accurate to resolve separate peaks for the two 
frequencies, so the peaks would merge to broaden the 
measured resonance width. However, the peak for uj + is 
likely to be higher than the peak for uj- owing to greater 
probability for ejection from the beam. Hence, the center 
of the merged distribution will be displaced to a higher 
frequency. If this explanation is correct, then an increase 
in experimental resolution will separate the two peaks, 
and their relative heights will measure the relative prob- 
ability of ejection at the two frequencies. 

Though the string approximation to the crystal po- 
tential (185) is useful for semiquantitative analysis of 
zitter resonance, as we have seen, it breaks down com- 
pletely at radii near the screening radius a = 0.190 A; 
for then, as Lindhard [26] has shown, the electron's colli- 
sion time with a nucleus is comparable to the travel time 
between atoms. In picturesque terms, the electron be- 
gins to "feel" individual atoms rather than a continuous 
string. Within this domain, our analysis of electron mo- 
tion remains qualitatively the same, but a more realistic 
crystal potential is needed for accurate quantitative esti- 
mates. A hard lower limit on the radius is determined by 
the mean radius of nuclear zeropoint vibrations, which 
Debye theory estimates as 0.05 A [25]. 



If the idea of zitter resonance is taken seriously, there 
are many opportunities for new theoretical and experi- 
mental investigations. Increasing the resolution by three 
orders of magnitude will open the door to refined studies 
of frequency shifts, line splitting, spin effects and Zee- 
man splitting, all of which are inherent in zitter theory 
[30]. As has been noted, the most straightforward pre- 
diction of the zitter model is a second order resonance 
near 161.7MeV/c. In a first approximation, it can be an- 
alyzed in much the same way as here, though removal of 
electrons from the center of the beam may be enhanced 
by such processes as pair creation. 

A more detailed analysis of zitter resonance in channel- 
ing requires close attention to experimental conditions, so 
that will be addressed elsewhere. Classical particle mod- 
els have long been used for channeling calculations with 
considerable success. Besides being simpler and more 
transparent than quantum mechanical models, they even 
often give better results at high energies. For present 
purposes, the zitter model differs from the usual classi- 
cal model only by the zitter dipole interaction. As we 
see in the next Section, though zitter is inherent in the 
Dirac electron theory, it is unlikely that the channeling 
resonance effect can be derived from the Dirac equation 
without some modification such as projection into a Ma- 
jorana state. 



F. Solution by separation of variables 

An obvious strategy for solving the rotor equation of 
motion is to separate internal zitter oscillations from the 
effect of external field on overall motion of the comoving 
frame. Spinor equations admit a natural way to do that. 
For a given electromagnetic field F, the rotor R = R{t) 
can be factored into 



R = LU, 

where rotor L satisfies the equation 

q 



L 



2m e 



-FL. 



(213) 



(214) 



Moreover, instead of using zitter averages as before, we 
can define a mean velocity v for the zitter center by re- 
quiring 



R^oR — LjoL so Uj U = '-fQ. 



(215) 



This generalizes (144) by dropping the requirement that 
L be a boost. However, the boost condition placed no 
restrictions on equations for the velocity, whereas here 
we have 



m e v = qF • v, 



(216) 



for which the relation to the momentum equation (95), or 
even its zitter mean (135), is problematic in general. This 
drawback aside, the approach has sufficient advantages to 
consider its use for approximate solutions at least. 



20 



Its first big advantage is evident in the transformed 
spin equation of motion (163), for which the fl has been 
adjusted to cancel F , so it reduces to 

5 , o = ^^-(ei+ie 3 )=me 2 . (217) 

Z UT 

This is equivalent to the rotor equation 

1 1 
U = -une 3 U = -Uut 3 u> 

with e fe = UcT k U. (218) 

It has the obvious solution U = Uq exp ier 3 (p, though w = 
if must be obtained from other considerations. 

A second big advantage is that exact solutions of (214) 
have already been found for several kinds of electromag- 
netic field [1, 31]. Unfortunately, they all presume that 
motion is along a timelike curve x(t), whereas, in the 
present case F = F(z) = F(x+r) is given along the light- 
like curve z(t), so problems arise in relating the curves 
x(t) and z(t). 

Another difficulty is that the constraint equation relat- 
ing zitter phase ip to proper time is nonlinearly related 
to the rotor equation of motion, as is obvious when ex- 
pressed in the form 

ip = uj e + — < FLSqL > . (219) 

m e 

Consequently, a closed solution of the coupled equations 
is not to be expected in most cases. However, the in- 
teraction term in (219) is invariably much smaller than 
the constant term oj e , so it can often be neglected or, in 
any case, incorporated by iteration with the equation of 
motion. 

The complications just described disappear when F is 
uniformly constant, as in the case of any combination 
of homogeneous electric and magnetic fields. Without 
working out the details, let's note some special features 
of this case that can contribute to an elegant solution. 

For constant F, equation (214) has the simple solution 

L.exp(^-Fr). (220) 

Since FL — LF in this case, (219) reduces to 

tp = lu e + — < FS >= oJ e + — (E • ei — B • e 3 ). (221) 

m e m e 

The initial phase can be chosen to make the angle depen- 
dence of the zitter oscillations explicit: 

ip = uj e - — B • e 3 + — IE X e 3 | sin<£. (222) 

This equation is akin to Kepler's equation in celestial 
mechanics relating incommensurate time variables. Since 
the oscillating term is so small, solution by perturbation 
expansion is quite satisfactory. 

Furthermore, from the momentum equation of motion 
(95) we find immediately that p — qF • z is a constant 



of motion. And the equation of motion (101) for tt = 
p — m e u reduces to 

7r = Q • tt = —mu) e ei, (223) 

so tt 2 = p 2 — 2m e m is another constant of motion. 

Finally, a nice trick for integration to get an explicit 
equation for the particle history is given in [31]. This is all 
we need for the complete solution of motion in arbitrary 
homogeneous fields. Note that complications of solving 
for the wobbling momentum vector have been avoided. 

VIII. ZITTER RESONANCE AND FIELD 
THEORY 

Like the original Dirac theory, our zitter model de- 
scribes a test particle, which is to say that the electro- 
magnetic field that it generates is ignored. But if we take 
the zitter charge circulation seriously, it is evident that 
the particle must be the source of a fluctuating electro- 
magnetic field. The average of that field over a zitter pe- 
riod will give us the static field of a magnetic dipole and 
an effective Coulomb field with a virtual source at the zit- 
ter center. Superimposed on these mean values we have 
an electric dipole field rapidly rotating with the zitter fre- 
quency. This picture raises many questions that cannot 
be answered without constructing a complete electromag- 
netic field theory with zitter sources. That is too much 
to attempt here, but we can identify some of the issues 
and prospective physical implications. 

We consider here only a classical field theory, leaving 
issues of second quantization and many particle theory 
for another day. The first question is: What maintains 
the zitter charge circulation? Classical electrodynamics 
implies that all accelerating charges radiate, which lead 
to collapse of circulation. An obvious suggestion is that 
circulation is maintained by self interaction of the charge 
with its own field! Of course, that cannot be regarded as a 
satisfactory answer until the self-interaction is calculated, 
and no such calculation has ever been successful in any 
version of classical theory. Nevertheless, it does provide 
a new constraint on the self-interaction problem, namely, 
that its solution should account for electron spin as well 
as mass. Besides that, there are tricky questions about 
singularities in the field of a lightlikc source. Evidently 
some new idea will be needed to make zitter field the- 
ory work. However, we can take an intermediate stance, 
presuming the zitter field as real and exploring its impli- 
cations. 

The dipole interaction d • E is, of course, unobserv- 
able for laboratory electric fields, because d rotates so 
rapidly. However, if E is the oscillating zitter field of 
an electron, the possibility of resonance exists. The fol- 
lowing speculative survey suggests that many surprising 
features of quantum mechanics might be consequences of 
zitter resonances. 

Electron diffraction as resonant momentum trans- 
fer mediated by the zitter field. An electron incident on 
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a crystal is preceded by its zitter field which is reflected 
back from a periodic array of scattering centers. At the 
Bragg incidence angles, phase coherence of the reflected 
field oscillations with the electron zitter provides a mech- 
anism for resonant momentum transfer. The propagating 
zitter field has a wavelength that may be identified with 
the electron de Broglie wavelength in quantum mechan- 
ics. 

Quantized atomic states might be explained by the 
same mechanism as diffraction, that is, as resonances of 
the electron zitter with its own field reflected off the nu- 
cleus. This raises questions about the topology of the 
zitter field that may be related to Berry phase and the 
Aharonov-Bohm effect. Each electron in a multi-electron 
atom will have its own zitter frequency, shifted by its in- 
teraction energy in the atom. Consequently, an atom or 
molecule will have a zitter field composed of the multiple 
frequencies of its constituents, and it will have not a sin- 
gle de Broglie wavelength but a spectrum of wavelengths. 
Perhaps this spectrum will show up in diffraction exper- 
iments. 

The Pauli Principle as resonant phase locking. In 
general, the zitter oscillations of different electrons will 
not be coherent, because of doppler shifts and frequency 
shifts due to local interactions. However, when two elec- 
trons are in the same quantum state the necessary condi- 
tions for resonant phase locking are fulfilled. It remains 
to be proved that antiparallel alignment of spins is nec- 
essary for a stable state. 

Van der Waals and Casimir forces are generally 
explained as effects of coherent coupling between fluctu- 
ating dipoles. This involves coherent coupling of phases, 
which is what zitter is all about. The zitter field serves 
as a mechanism for maintaining that coupling that may 
have additional consequences. As the dipoles are mutu- 
ally induced, the interaction mechanism appears to be 
essentially the same as in electron diffraction. 

The Lamb Shift is commonly attributed to a smear- 
ing out of electron position due to vacuum fluctuations. 
Obviously, the zitter offers an alternative explanation for 
smearing out. Indeed, one can conceive an electromag- 
netic vacuum field composed of a stochastic combination 
of zitter fields of all existing charged fermions. 

Similar speculative explanations can be adduced for 
tunneling, anomalous magnetic moment, covalent bond- 
ing and other important quantum effects. Even Bose- 
Einstein condensates might be explained by resonant zit- 
ter coupling of constituent fermions. Zitter field theory 
must provide reasonable accounts of all these phenom- 
ena before it can be regarded as a viable extension of 
standard quantum mechanics. However, we don't have 
to wait for that! If zitter is a real physical phenomenon 
we should be able to devise experiments to observe it di- 
rectly. In fact, as we have noted, that may already have 
been done! 



IX. ZITTERBEWEGUNG IN DIRAC THEORY 

This section describes the intimate relation of the zitter 
model to the Dirac equation. The relation should not be 
surprising, as it was study of Dirac theory that led to 
the zitter model in the first place. Conversely, we shall 
see that the zitter model suggests modification of Dirac 
theory to incorporate deeper zitter substructure, with 
new physical implications including a surprising possible 
connection to weak interactions. 

Our first task is to match up variables and dynami- 
cal equations in the zitter model with observables and 
dynamics in Dirac theory. An exact match is not to be 
expected, as the particle based zitter model can only be 
related to the Dirac field theory by some sort of averag- 
ing or projection process. Nevertheless, the comparison 
reveals specific similarities and differences that must be 
addressed in establishing a firm connection between zit- 
ter model and Dirac theory. This leads to suggestions 
for modifying the Dirac theory and further research to 
resolve outstanding issues. 

In the language of STA, the real Dirac equation has 
the form 

7^(^7271^-9-4/^) = m oi>lo, 

or \7tpier 3 h ~ qAip — m e ip-f > (224) 

where the Dirac wave function is a real spinor field 

$ = ^(x) = {pe lP ) 1/2 R. (225) 

This version of the Dirac equation is fully equivalent to 
the standard matrix version [1] , but it has great advan- 
tages for analyzing the structure of the Dirac theory as 
shown in the following. 

The Dirac wave function determines a frame field of 
local observables 

= pe^, where e p = R-y^R = e M (x) (226) 

and p — p(x) is interpreted as a scalar probability density, 
in accordance with the interpretation of the Dirac current 
ipjoip — peo as a probability current. The vector fields 

v = e = RjqR = v(x), (227) 

and 

s = ie 3 = \hR^R = s(x) (228) 

are interpreted as local velocity and spin observables for 
the electron. Note that these quantities are algebraically 
identical to the expressions for velocity and spin in the 
zitter model. Likewise for the other observables 

u = e + e 2 = R-/+R = u(x), (229) 

and 

S = isv = \Ricr s hR = S(x). (230) 

The difference is that the observables here are vector and 
bivector fields, whereas in the zitter model they are de- 
fined on a particle history. Our next task is to compare 
them dynamically. 
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A. Zitterbewegung along Dirac histories 

The mass term in the Dirac equation (224) can be writ- 
ten in the form 

m e ~ip~f = m e (pe ip ) 1/2 R~/ RR = m e ve~ if3 i/j. 

Whence, the Dirac equation can be reformulated in the 
compact form 

Dip = j"D^ = 0, (231) 
where a gauge invariant coderivative is defined by 

£>„ = d„ + £fi„, (232) 
with the bivector-valucd connexion 

fl^ = ^{m e v^e- lfi + qA^ei. (233) 

The purpose of introducing this coderivative is not to re- 
duce the Dirac equation to the maximally compact form 
(231), but to reveal that all the essential physics is con- 
tained in the connexion M . 

Our aim is to compare zitter dynamics along the 
history of a zitter center with Dirac dynamics along 
"streamlines" of the Dirac current. To that end, we eval- 
uate the directional coderivative along a Dirac history as 
follows. From 

£>V>7o = D{ve- 1 ^) = \D(vip) - i(V/3)W]e~ l/3 = 

we obtain 

D(vtP) =2vDil) + (Dv)tP = i(V(3)vip. (234) 

This separates into the familiar conservation law for the 
Dirac current and a dynamical equation for rotations 
along a Dirac streamline: 

p + pD • v = D • (pv) = V • (pv) = 0, (235) 

v ■ DR = R+ ±Q(v)R = -\[D A v + i(v A Vf3)]R. (236) 

The overdot indicates the directional derivative v • V and 
(233) gives us 

Q(v) = (u cos/3 + ^j-A • v)e 2 e 1 + (w sin/3)e e 3 (237) 

D Av = V Ad| (w sin/3)e 3 e . (238) 

When these are inserted into (236), the boosts in the eoe^ 
plane cancel to give us 

R = ^[(u) cosp+^-A-v)e 1 e 2 -VAv-i(vA\'(3)}R.(239) 

This is an exact result. It does indeed exhibit the familiar 
zitter rotation in the e 2 ei-plane, though the frequency 
seems to be different. 

We still need to evaluate the curl of the velocity field to 
appreciate its effect on the dynamics. A general expres- 
sion has been derived elsewhere [2, 3], but our purpose 



here is served by the eikonal approximation, expressed 

by 

VV>727i7i = V(^oe~ 7271 * /?i 727i?i) = (V$)V- (240) 

Inserting this into the Dirac equation (224), we obtain 

V$ = m e ve- lP - qA. (241) 

This implies e 1 ^ = ±1, where the choice of sign depends 
on the chosen sign of charge. We adopt that approxi- 
mation only in this equation, as the parameter (3 is too 
important to ignore completely. Then the curl of (241) 
gives us 

-m e V A v = qV A A = qF. (242) 

This is equivalent to the Lorentz force equation for a fluid 
of charge with uniform density, as seen by "dotting" with 
v, to get 

v • (V A v) = v • Vu = — — v ■ F = —F ■ v. (243) 

More generally, we insert (242) into (236) and use (238) 
to get the spinor equation of motion 

R= ±[uj v e 1 e 2 + — F-t(vAVP)]R, (244) 

with lo v = loq cos [3 + (2q/h)A • v. This equation must be 
compared with the analogous expression (138) for rota- 
tional velocity in the minimal model. The general form 
is very similar, but we are not equipped to account for 
the differences for reasons to be discussed. 

The zitter frequencies in the two equations are not 
equivalent, but they do have the same free particle limit 
u>o- The bivector i(v A V/3) also contributes to the rota- 
tion rate, as implied by 

v • [i(v A V/?)] = -i(v Av A V/3) = 0. 

However, the physical significance of this term remains 
obscure. The apparent absence of a Stern-Gerlach force 
in the Dirac version (244) is noteworthy, but we cannot 
be sure that it is not buried in terms that we do not 
understand. 

It might be thought that the Dirac equation is more 
fundamental than the zitter model because interaction 
comes from the vector potential A alone and interaction 
with the field F arises only indirectly, for example in the 
manner described above. A famous consequence of this 
is the derivation of the gyromagnetic ratio g = 2. In con- 
trast, the Lagrangian (59) for the zitter model appears 
to presume the electron magnetic moment with separate 
coupling constants for A and F interactions. Note, how- 
ever, that precisely two independent constants are pre- 
sumed in both models. The rest mass is presumed in the 
Dirac equation, but that is replaced by the zitter radius 
in the zitter model. It remains to be seen which is more 
fundamental. 
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B. Zitterbewegung substructure 

In preceding sections we saw that a rotating electric 
dipole is the hallmark of zitterbewegung, so one wonders 
why it has attracted so little attention in accounts of 
Dirac theory. In the original paper introducing his equa- 
tion [32], Dirac concluded that the electron has both a 
magnetic and an electric moment, the magnetic moment 
being the same as in the Pauli model. However, he said, 
"The electric moment, being a pure imaginary, we should 
not expect to appear in the model. It is doubtful whether 
the electric moment has any physical meaning." 

It is worth translating Dirac's argument into STA, es- 
pecially since its mathematical content has been retained 
in the current literature. One simply "squares" the op- 
erator on the wave function in (224) to get 

{-h 2 V 2 + q 2 A 2 )tjj - 2q(A ■ V + %F)ilntr 3 h = m 2 e tp ,(245) 

where F = VA = V A A with the Lorenz condition 
V • A = 0. This is the Klein-Gordon equation with an 
extra term that explicitly shows the action of bivector 
F "rotating" the wave function. The interaction energy 
density associated with this term is proportional to 

\ < F^ia 3 h^ >=< FSpe i0 > 

= -p(B v • s cos/3 + E^ • s sin/3), (246) 

where the u-split introduced in (142) has been used on 
the right hand side. The B„ -s cos /3 term is recognized as 
the Pauli term except for the strange cos/3 factor. The 
E„ • s sin /3 term is what Dirac identified as an imagi- 
nary dipole moment. He never mentioned the electric 
dipole again. In his influential textbook [33] he simply 
suppressed the offending term by a subterfuge advertised 
as a change in representation. Then he killed the term 
with an approximation that amounts to sin/3 = and 
never looked back. We find an alternative resolution of 
this "dipole problem" below. 

Dirac was soon convinced by Schrocdingcr [5] that zit- 
terbewegung is foundational to electron theory and he 
argued the case vigorously in his textbook [33]. As his 
argument is still widely accepted [34], it deserves com- 
ment here. 

Dirac introduces a position operator by identifying a 
velocity operator as its time derivative, and he followed 
Schrocdinger in integrating the equation for the free par- 
ticle case. He identifies his at matrices as velocity op- 
erators and claims that their eigenvalues ±1 correspond 
to measured values of electron velocity, asserting: "we 
can conclude that a measurement of a component of the 
velocity of a free electron is certain to lead to the re- 
sult ±c." (Dirac's italics) From the STA point of view 
this argument and its implications are bogus, for reasons 
explained elsewhere [1]. However, we agree with the as- 
sumption that the electron moves with the speed of light. 
We differ in identifying the local electron velocity with 
the null vector u = R^ + R introduced above. 



We also agree with Dirac in attributing the origin of 
spin to zitter. Dirac concludes: " Our argument is valid 
only provided the position of the particle is an observ- 
able. If this assumption holds, the particle must have a 
spin angular momentum of half a quantum." However, 
Dirac's analysis of zitter and spin never went beyond the 
free particle case. He overlooked (or dismissed) the infer- 
ence that his account of spin arising from charged parti- 
cle circulation implies a rapidly rotating electric dipole. 
Though zitter obviously arose from wave function phase 
oscillations in Schroedinger's free particle analysis, Dirac 
never considered a general connection of zitter circulation 
to wave function phase. The present study can be re- 
garded as an extension of Dirac's analysis to incorporate 
these features in a general theory of zitter in quantum 
mechanics. 

As the Dirac equation has an unsurpassed record of 
success in QM and QED applications, it is imperative to 
reconcile it with any proposals about zitter. If the zitter 
model describes substructure in electron motion that is 
not captured by the Dirac equation, it must at least be 
related to the Dirac wave function by some sort of aver- 
aging process. Without attempting a definitive reconcil- 
iation, let us note some issues that must be addressed. 

We know that the conserved Dirac current pv deter- 
mines a congruence of curves (or streamlines) for every 
solution of the Dirac equation. As Bohm and Hiley have 
argued at length [4], each of these curves x — x(t) can 
be regarded as a possible path for the electron weighted 
by a relative probability p = p(x) that the electron ac- 
tually followed that path. This is a viable particle inter- 
pretation of quantum mechanics. However, a refinement 
is necessary to account for zitterbewegung, which sug- 
gests that the actual particle paths are lightlike helices 
with tangent vector u = u(x) at each spacetime point. 
The simplest refinement would have each of these light- 
like paths winding around a Dirac streamline, but this 
possibility is questionable without deriving it from the 
zitter model, at least approximately, by a well defined 
averaging process (yet to be determined) ! 

A crucial problem is to justify the weighting of paths 
by the probability density p = p(x). A new possibility 
is suggested by the fact that the analogous quantity in 
our zitter particle model is a timescale factor p = p(ip) = 
J^, which, as we have noted, would arise naturally from 
averaging over paths with different proper times but a 
common phase angle tp. Thus, the putative probability 
density in the Dirac equation might be derivable as a time 
scale weighting on a congruence of particle paths! In any 
case, time scaling in the zitter model must be reconciled 
with the probability interpretation in the Dirac theory. 

It seems likely that a suitable averaging process relat- 
ing particle histories to the Dirac equation will involve 
time averaging as well as ensemble averaging with con- 
straints. Analogy with the zitter model already suggests 
that Dirac observables for velocity v and spin S corre- 
spond to zitter time averages, but the phase of the Dirac 
wave function is directly comparable to the zitter phase. 
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The Schroedingcr equation is a nonrelativistic approx- 
imation to the Dirac equation that freezes spin but pre- 
serves zitter oscillations in the phase of the wave function 
and in coupling of phase to amplitude, as specified, sur- 
prisingly, by the mysterious parameter (i [31, 35]. As 
shown in these references, assuming (3 — completely 
decouples phase from amplitude and so eliminates all 
QM effects. Clearly, therefore, the role of (3 must be 
accounted for in any averaging process. 

The fact that the kinematic state of a particle with 
zitter is described by a rotor while the Dirac wavefunc- 
tion is also a spinor suggests that the QM superposition 
principle can be construed as an average over particle 
rotor states. The undeniable success of the Schroedingcr 
equation suggests that, approximately at least, the super- 
position is an average over phase factors, such as Feyn- 
man's sum over paths. This certainly produces the cou- 
pling of phase to amplitude so characteristic of QM. Note 
that an average over rotors blurs any average over paths, 
because path velocity is a bilinear function of the path 
rotor. Finally, it should be mentioned that gauge invari- 
ance provides a strong constraint on assembly of phases 
for distinct particle paths into a coherent ensemble. 

Clearly, deriving the Dirac equation from zitter sub- 
structure is a nontrivial problem. Happily, we don't have 
to wait for a solution to make progress in studying the zit- 
terbewegung. Here is a promising alternative approach: 

C. Putting zitterbewegung into the Dirac equation 

We have seen that physical interpretation of the Dirac 
equation is crucially dependent on identification of a par- 
ticle velocity observable, which requires theoretical as- 
sumptions beyond the Dirac equation itself. Historically, 
the mass term in the Dirac equation led to the conserva- 
tion law for the Dirac current and its interpretation as a 
probability current, with an implicit identification of par- 
ticle velocity. In the STA version of the Dirac equation 
(224), explicit appearance of the vector 70 shows that 
the velocity vector v = RjoR is inherent in the structure 
of the equation. However, we have identified the vector 
u = R / y + R as a better candidate for electron velocity. 
This suggests a slight modification of the Dirac equation 
to replace v by u. Accordingly, we change the mass term 
m e V>7o to 

TO e -0|(7o + 72) = m e -0+7o, (247) 

where 

^+=^±(1 + <t 2 ). (248) 

In common parlance, this is a projection of a 4- 
componcnt spinor -0 into a 2-component spinor ip + . As 
only these components are now relevant to the electron 
velocity, we should perform the same projection on the 
other components of the Dirac equation. That requires 
modification of the vector potential term to achieve a 



2-component equation. Thus we arrive at the modified 
Dirac equation 

Vip + icr 3 h - qAip + (T 3 = TO e ^+7o • (249) 

Let's refer to this as the zitter Dirac equation. Note that 
it is invariant under projection from the right by i(l — 

To be assured that this modification of the Dirac equa- 
tion has not damaged its essential physical meaning, we 
note that the zitter free particle solution (57) is also a 
solution of the Dirac equation, which, when substituted 
into (249) with arbitrary initial conditions and i = ier 3 , 
yields the algebraic relation 

pu = m e (l — e2eo). (250) 

This is identical to the relation between momentum and 
velocity found in (80) for the zitter model. In contrast 
to the prosaic relation p = m e eo from the free particle 
solution to the ordinary Dirac equation, this relation in- 
cludes zitter in the vector e 2 , as it rotates with the zitter 
phase. Thus, the form of the zitter Dirac equation brings 
the null velocity observable to the fore. 

It is also readily shown that the zitter Dirac equation 
has the same electromagnetic gauge invariance as the or- 
dinary Dirac equation, though we shall see that the gauge 
generator assumes a different form when the gauge group 
is generalized to incorporate electroweak interactions. 

Now check the observables. With respect to the pro- 
jected wave function, the observables (229), (230) become 

^+7+^+ = 2V>+7o^+ = V'T+V' = PU, (251) 

and 

ip + icr 3 frip + = ip±(l + a 2 )i<T 3 hi! 

= ^7+7iV> = p\ue-^e x = pS. (252) 

Thus we get the same zitter velocity u, but the spin bivec- 
tor S is replaced by a null spin bivector S, which is iden- 
tical in form to the spin (73) in the point particle model, 
and, as in that case, the duality factor can be absorbed 
into a rotation of the vector e\. We see immediately that 
the interaction energy density (246) becomes 

I < Ftp + icr 3 fi^ + >=< FpS > 

= p(E v • d - B„ • s), (253) 

in perfect accord with equation (143) of the zitter model. 
Thus, the change in observables by projection on the 
wave function appears to eliminate the parameter (3 and 
its problems of physical interpretation. 

Further insight comes from the following Lagrangian 
for the zitter Dirac equation: 

LzD = 2 ([-KVip+i'Yo 

+qAip + ~/ + m e ip + (T 3 ]^) . (254) 
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Note that the interaction term has the usual form A • 
J =< A J >, but now the charge current is a null vector 
field 

J = 2#+7o^ = #7+V> = 1P U , (255) 

as expected. It follows from the zitter Dirac equation 
that this current is conserved, though it is not conserved 
in standard Dirac theory. 

The usual Dirac current is not obtainable as a bilinear 
observable of the wave function ip + . However, as noted 
before, it can be obtained as a zitter average v = u of the 
zitter velocity. This suggests that one should try to derive 
the zitter Dirac equation, rather than the Dirac equation 
itself, from the zitter model. But that possibility will not 
be explored here. 

The zitter Dirac equation offers a new perspective on 
the significance of negative energy in Dirac theory. Re- 
call that negative energy solutions were first regarded 
as a serious defect of the Dirac equation. Schrocdingcr 
showed that they are essential for a general solution even 
in the free particle case when he "discovered zitterbe- 
wegung" as interference between positive and negative 
components of a wave packet. In an audacious effort 
to save the theory, Dirac identified the negative energy 
states with an "anti-electron" and invoked the Pauli prin- 
ciple to suppress them (hole theory). Miraculously, the 
positron was discovered shortly thereafter, so the defect 
was transmuted to a spectacular triumph! 

A standard conclusion from all this is expressed by the 
following quotations [34]: "The zitterbewegung demon- 
strates in a real sense a single particle theory is not possi- 
ble." "The difficulties with the negative energy states of 
the Dirac equation almost of necessity demand a many- 
body theory." "Hole theory is a many-body theory de- 
scribing particles with positive and negative charge. The 
simple probability interpretations of the wave functions 
acclaimed in a single-particle theory cannot be true any 
longer, because the creation and annihilation of electron- 
positron pairs must be taken in account in the wave func- 
tion." In other words, quantum field theory is needed to 
explain zitterbewegung! 

For a different perspective on negative energy we note 
that tp c = iper 2 is the charge conjugate solution of the 
real Dirac equation [1]. Hence we can cast the zitter wave 
function (248) in the form 

^i(^ + ^ c ). (256) 

This expresses Dirac's negative energy solution as an es- 
sential component of the zitter rather than an antiparti- 
cle. It is an alternative splitting of Dirac's 4-componcnt 
wave function into a pair of 2-component wave functions 
for different particle states. The physical issue is this: 
Which components of the Dirac wave function should be 
identified with the electron? The zitter component ?/>+ 
describes an electron with zitter motion. Standard quan- 
tum field theory splits the zitter into positive and neg- 
ative energy components and then reassembles it later 
from pair creation and annihilation. In most calculations 



the end result will be the same, because both approaches 
start from the same Dirac equation. Contrary to the 
standard Dirac equation (224), the zitter Dirac equation 
(249) is consistent with a single particle model of zit- 
ter without the strenuous expedient of field quantization. 
This is not to deny that some version of quantum field 
theory is necessary to account for creation and annihila- 
tion of particles. The problem is to devise experiments 
that identify the basic particle states. 

D. Electroweak interactions 

We have seen that only half the Dirac wave function 
■0+ = tp^(l+eT 2 ) is needed to describe the electron. What 
can be said about the other half ip- = — o^)? An 

attractive answer is suggested by gauge theory. 

In the real Dirac equation (224) an electromagnetic 
gauge transformation of the wave function is multiplica- 
tive on the right with the form 

^ -» V" = ipu, (257) 

where U — exp(z<X3x). We look to adapt this transfor- 
mation to the zitter Dirac equation (249) and generalize 
it in a way that preserves essential structure of the equa- 
tion. To preserve the mass term, we require 

UjoU = 7o- (258) 

The general solution of this equation has the form 

U = e* ie e* ix , (259) 

where 9 = 8\<J\ + 82&1 + 83(73. Remarkably, this is 
the gauge group SU(2)<8>U(1) of electroweak theory. It 
strongly suggests that the geometric structure of elec- 
troweak theory is already inherent in Dirac theory! In- 
deed, it requires that we identify = ip v with the neu- 
trino, just as we identify ip + = ip e with the electron. 

Thus, we construe the Dirac wave function as a lepton 
wave function composed of electron and neutrino compo- 
nents: 

V = ^(l + <r 2 ) + V4(l-<r 2 )=^ e + Vv (260) 

To adapt the Dirac equation to this interpretation, we fol- 
low standard electroweak theory in introducing a gauge 
covariant derivative 

with Wfj, = ^(ffA^ - g'B^). (261) 

In the standard way, gauge covariance is assured by re- 
quiring that the electroweak connexion satisfies the trans- 
formation law 

Wn - W; = UWpU - Ud^U. (262) 

Accordingly, the Dirac equation generalizes to the lepton 
wave equation 

^V^icr^h = TOeV'eTo- (263) 
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For the electron component, this reduces to the zitter 
Dirac equation (249) when weak interactions are turned 
off. 

Details of this model integrating zitterbewegung with 
electroweak theory are discussed elsewhere [30]. It is 
amusing to recall that one of the early suggestions to 
account for extra components in the Dirac wave function 
was to identify them as proton states. Long dismissed 
because the particles had different mass, the idea returns 
again as an electron-neutrino wave function. 

X. CONCLUSIONS 

The zitterbewegung, if it turns out to be physically 
real, is belated confirmation of de Broglic's original hy- 
pothesis [36] that the electron has an internal clock with 
period precisely equal to twice the zitter period, precisely 
the relation between the period of a rotor and that of a 
vector it rotates. 

As we have seen, the physical signature of zitter is a 
rotating electric dipolc with ultra high frequency If this 
exists, its implications for quantum mechanics will be 
far-reaching. Evidently it can be incorporated in Dirac 



theory by subtle changes in the specification of observ- 
ables and the structure of the Dirac equation. 

Experimental confirmation of the zitter should stim- 
ulate research on its proposed incorporation into elec- 
troweak theory. Then study of zitter self-interaction 
should look for excitations explaining the three lepton 
families. Finally the strong analogy between electroweak 
interactions of leptons and quarks suggests that one 
should investigate modifications of zitter structure to 
model quarks and strong interactions. All this should 
go hand-in-hand with development of zitter field theory 
and reconciliation of it with quantum field theory. 
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